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DYNAMIC ANALYSIS AND RESPONSE OF 
AIRCRAFT ARRESTING SYSTEMS 


Robert S. Ayre, A.M. ASCE 
and Joel I. Abrams,2 J.M. ASCE 
(Proc. Paper 1580) 


SYNOPSIS 


A general analysis of aircraft arresting systems, treating them as lumped 
parameter systems and omitting consideration of the effects of wave travel 
in the cable, is presented. The differential equations for aircraft displace- 
ment and for cable elongation are solved in a step-wise manner by phase- 
plane graphical methods. 

The systems are highly non-linear owing to geometry, to large displace- 
ments, and to the characteristics of the energy dissipators. Special linear 
cases have been shown for comparison. The investigation includes studies 
of the effect of varying the frontal-shape and rise-time of the energy 
dissipation function, the initial cable tension, the mechanical advantage in 
the cable and sheave system at the arresting engine, the deck pendant length, 
and the mass of the moving parts in the arresting engine. The results are 
also of interest in the general problem of arresting the motion of a moving 
mass. 

Numerical examples and a brief comparison with experimental results 
are included. Response spectra of maximum cable tension are presented 
for a wide range of operating conditions. 


INTRODUCTION 


The purpose of this paper is three-fold, first to present a “lumped para- 
meter” analysis of the dynamics of aircraft arresting systems, second to 
show applications of the phase-plane step-wise graphical method of solution 
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to the transient dynamics of highly non-linear mechanical systems, and third, 
to show the spectra of maximum response for a family of arresting systems. 

Aircraft arresting “gear” constitute an interesting class of dynamical 
systems that include many fundamental problems in mechanics. They consist 
primarily of cables and pulleys for engaging the landing airplane and trans- 
mitting the arresting forces to an energy dissipator, the arresting engine. 
Since the airplane is brought to rest in about one second, the rate of energy 
dissipation is high and the inertia forces in the system are large. Arresting 
systems are used not only on board aircraft carriers but also sometimes on 
land. The mechanical features of the systems may differ considerably. 

The motion of the airplane during arrestation is essentially non- 
oscillatory and may be determined to a good approximation by a simple 
energy solution. The cable tension, however, consists of oscillatory com- 
ponents as well as of the “static” component, the oscillations arising be- 
cause of the inertia force concentrated at the arresting engine and also be- 
cause of the distribution of mass along the cable. 

This analysis is based on the differential equations of motion of an 
assumed “lumped parameter” system, and omits consideration of the effect 
of elastic wave travel in the cable. The latter has been investigated by 
F. O. Ringleb(1,2,14) and others,(3) and currently is being investigated by the 
writers using a small scale mechanical mode}l.(15, 16) 


Basic System 


An elementary arresting system, shown in Fig. 1, consists of a cable and 
sheave system and arresting engine for dissipating the kinetic energy of the 
aircraft. The analysis is limited to a symmetrical system. The arresting ie 
engine, represented by the dissipation function F, may be a hydraulic- ie 


cylinder-type energy dissipator or a dry friction brake. The symbols are 
as follows: 


x displacement (runout) of aircraft, 
R maximum displacement (maximum runout) of aircraft, 


x,X aircraft velocity and acceleration, 


Vo velocity of aircraft at impact with cable (initial velocity), 


s,§ displacement and velocity of energy absorption element in arrest- 
ing engine, 


Ss maximum value of s, 

t time, 
M = . mass of aircraft (assumed to act as a rigid body), 4 
ni = = mass of concentration at energy absorption element, a 
u mass of cable per unit length, = 
A effective cross-sectional area of cable, 

E effective modulus of elasticity of cable, 


ton 
Sas 
A 
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one-half length of cable, 
€ elongation of one-half of cable length, 
L one-half length of “deck pendant” cable, 
F dissipation force at arresting engine; in general, a function of t, 
s, 8, and S raised to some power, 
N “mechanical advantage” in pulley system at arresting engine. 
“cable runout,” 
cable tension, 
To initial cable tension. 


I. Elasticity and Distributed Mass of Cable Neglected 


Energy Solution 


An energy solution results in a good approximation for aircraft displace - 
ment because the displacement is mainly a function of the initial kinetic 
energy, of the geometry of the system, and of the dissipation. The energy 
solution is useful not only as a basis for overall design but also as a first 
approximation for the dynamical solution. 

If secondary fructional losses and the storage of elastic energy are 
neglected, the basic energy relationship is given by 


Ss 9 ? 
flr as =4 - we? (1a) 
0 
consequently, 
Ss 
f F ds = 5 WV (1b) 


Since the elongation of the cable is assumed to be negligible, the relation- 
ship between the aircraft displacement and the arresting engine displacement 
is as follows: 


1 + L- (2a) 
consequently, 

ret (2b) 

N ( VR - L) 


Differentiating s with respect to x, one finds 


from which the following relationship between engine velocity, aircraft dis- 
placement and aircraft velocity may be obtained: 


ds x 
dx 
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The differential equation of motion of the airplane is given by 


Mx = - 27 —*~—— (4) 


Constant dissipation force. A “rectangular” force-displacement relation- 


ship (Fig. 2) results in the smallest maximum value of the non-oscillatory 
component of cable tension. Equation (1b) then becomes 


and by applying equations (2b) and (3) the following expression for cable ten- 
sion is found: 


c  2N (6) 
ac Vee +L" - 


Curves of T,./W versus R/L have been shown in Fig. 3. 
From equation (la) it can be shown that the velocity for this case is given 
by 


leo (7) 


+ 


e and from equation (4) that the acceleration is 


V x 


2 
Vr" +L” -L) 


In Fig. 4 the velocity and acceleration have been plotted as functions of the 
aircraft runout x for various values of R/L. 


Phase-plane Graphical Method of Solution 


Since equation (8) is non-linear in x, its integration, to result in an ex- 
pression for x as a function of time, presents some difficulty. The integra- 
tion may be carried out conveniently by a graphical step-by-step procedure 
in the phase plane”, as follows: 


*The phase-plane method has been fully described elsewhere; see References 
(5,6,7,8).. A step-by-step numerical method based on the equation of instan- 
taneous energies (eq. la) may of course be used. 


ae 
ae 
¥ 
. 
) 
Wee (2c) 
Ee Furthermore, the relation between cable tension and engine force is 7 
(3) 
\ [2 2 
vx 
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FIG.2 CONSTANT DISSIPATION FORCE 
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The general second-order differential equation may be written in the 
“6 form,” as shown by L. S. Jacobsen,(7) 


ef 2 
x +p «& -8) =0 (9a) 


where p is a constant of arbitrary value having the dimensions of “circular a 
frequency” (time-1), and where 6 is in general a function of x, x and t. It 
will be noted that the homogeneous equation, 


x +p x =0 (9b) 


is the differential equation for free vibration of an undamped linear oscillator. 
The displacement-velocity trajectory relating to equation (9b) is a circular 
arc with the center at the origin of coordinates (Fig. 5a). In the case of equa- 
tion (9a) the instantaneous trajectory is a circular arc, subtended by an angle 
p-dt, having its center on the x-axis at a distance from the origin equal to the 
instantaneous value of 6 (see Fig. 5b). 

If, instead of the infinitesimal dt, finite intervals of time At are employed, 
and if 6 is treated as a step-wise constant, the phase-trajectory becomes a 
succession of connected circular arcs of finite length, as in Fig. 5c, and the 
trajectory can be constructed with ordinary drafting instruments. This 
method of solution amounts to a step-wise linearization of the problem. 

The phase-trajectory determines velocity as well as displacement, and 
time may be determined by summation of the angles p-At. The construction 
of displacement-time and velocity-time curves may be made by direct projec- 
tion from the phase-plane. 


is not a simple function of engine displace ment and is given graphically. 
Writing the differential equation (4) in the form 


General dissipation force. In the following example the dissipation force 


x O (10a) 


F l (10b) 


Assume, for example, that F is a function of s as shown in Fig. 6a. By 
using equation (2a), F may be determined as a function of x, Fig. 6b, and 
assigning a value* to p, 6 may be plotted as a function of x, Fig. 7a. The 
phase-plane solution has been shown in Fig. 7b. The projections from the 
phase -plane, of displacement versus time and velocity versus time, have 
been shown in Figs. 8a and 8b. The curve of cable tension versus time ap- 
pears in Fig. 8c. Curves of engine displacement and engine velocity versus ‘ 
time may be plotted by use of equations (2a) and (2c). 


*In this example, we use p = 2V,/R. While p may be assigned any real 
value, it is advantageous to select it so as to result in desirable proportions 
of the phase-trajectory. 


‘le 

= | 
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Vx 
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X+p*x =O 


ja 


x+p*x = p*8 


= pAtnee 
= pAtne: 
pAtn 


FIG.5 ELEMENTARY RELATIONSHIPS IN 
THE 4 ,x PHASE- PLANE 
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The graphical construction can be made as accurate as one desires by em- 
ploying smaller time intervals. Furthermore, the time intervals need not be 
uniform. They may be varied in order to obtain greater or less detail in 
particular parts of the trajectory. 


Constant dissipation force. Comparing equation (8) with equation (9a), it 
is seen that for the case of constant dissipation force, 6 is given by the 
following expression: 


S=x . _1__ (11) 


2p + x" +L* 


The graphical solution follows the same procedure as for the case of the 
general dissipation force. 


Dissipation force as a function of engine velocity. If it is assumed that the 
dissipation force is proportional to engine velocity or to some power of 


engine velocity, the dissipation function takes the form 


F = cs)" (12a) 


where the constant C, must be chosen so as to satisfy equation (1b), that is 


Cc (12b) 
2fs (3) "ds 


(A preliminary value of the integral Sf (8)"ds may be determined by com- 


puting the area under an estimated curve of (8)" versus s.) An expression 
for 6 is found by substituting the above relationships in equation (10b) and 
by making use of equation (2c). 

For example, if n = 1, 6 is given by the following: 


5 (13) 
= x - 
+1) 


Only the end values of 6 are known, namely, at x = 0 and k = Vo, 6 =0 and 
at x = R andx =0, 6 =R. Intermediate values of 6 must be calculated as . 
construction of the phase-trajectory progresses. If the determined value of 

R differs from the required value, the difference may be corrected by revis- 

ing the first estimate of Cy, and repeating the phase-plane construction. 


II. Elasticity of Cable Included; 
Distributed Mass of Cable Neglected 


Two-Degree-of-Freedom System 


In the preceding section the elongation of the cable was assumed to be 
negligible. In this section we take into account the elasticity of the cable, 
but the mass of the cable is assumed to be concentrated at the point of contact 
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with the aircraft and at the arresting engine. Under these assumptions the 
aircraft and the arresting system constitute a two-degree-of-freedom system 
(Fig. 1) governed by the following differential equations of motion: 


for the aircraft: 


14 
(14) 
for the energy absorption element in the arresting engine: 

ms = - F + 2NT (15) 


Equation (14) is the same as equation (4). The cable tension T is a function 
of x and of s as follows: 


where € is the total elongation in one-half of the symmetrical cable system. 
Equations (14) and (15) must be solved simultaneously. This may be done by 
step-by-step numerical methods but they are very laborious. 


Reduction to a Single-Degree-of-Freedom System 


The aircraft displacement is largely independent of the cable elasticity 
and may be calculated quite accurately by simple methods of analysis. This 
suggests a means of simplifying the problem considerably without at the 
same time destroying its practical significance. The method is to calculate 
x as a function of time, neglecting the elasticity of the cable and using the 
methods of the previous section, and then to assume, in the solution of equa- 
tions (14) and (15), that x is known. This has the effect of reducing the sys- 
tem to one having a single-degree-of-freedom. 


Cable Elongation (or Cable Tension) as the Dependent Variable 


Beyond the elementary relationships determined by the energy solution, 
the primary concern is with cable tension rather than with a very accurate 
determination of the engine displacement. However, looking at equation (16) 
we see that the determination of cable tension may involve the calculation of 
small differences, and this is apt to require a large number of significant 
figures. Therefore, we look for a way of calculating the cable elongation ¢€ 
directly. In order to do this we will rewrite the problem in terms of € as 
the dependent variable. 

From equation (16) may be obtained the relation 


(17a) 


(17b) 


ae 4 

= 

x 

j 
a 

~ 

aA 

4 
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4 

= 

N L L 

and by differentiation with respect to time, 
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9 


The next step is to find expressions for x and x in terms of x and €. 
Making use of the differential equation of motion of the aircraft, equation (14), 
we find 


A relation between x and x may be found from equation (18) in the follow- 


so that x = yew . Then, from equations (16) and 


dx 
ing manner. Let v = at? = 


(18) may be obtained 


from which 


x? = SRE +k (VK (19) 


Differential Equation for Cable Elongation. By substituting expressions 
(18) and (19) into (17c) we obtain 


2 2 
9 a — 
+ 1°)° 
2 
_ 2ZAE €x _2 20 
(20) 


Rewriting equation (15), substituting T = AE € / gq » one finds 


and 

9 9 9 9 

9 9 
7 
Integrating once, we find = 
v dv - 1) dx 

x =v = (Ve - 

and substituting expression (17a) for s, we obtain 
a 
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_ F 2NAEE 
(21) 


Finally, by equating expressions (20) and (21) for ‘S, the differential equation 


for cable elongation is obtained. It may be written in the “ § form* for phase- 
plane solution, as follows: 


E+p'e = p29 (22a) 
where 
p- = 2NAE 
Ln (22b) 
and 
<n 9 5 (22c) 
P Ve +i" (x* + 


9 9 


2 
NM 


Frequency of Oscillation in Cable Tension. The constant p is the 
“circular natural frequency” of a single-degree-of-freedom system com- 
posed of the effective mass concentration m at the arresting engine and the 


2 
longitudinal spring constant k of the cable, where k = “ss 


The natural frequency f in cycles per second is given by 


(23) 
2r 27 V2m 
If the assumptions are valid we should expect to find, in a recording of 
cable tension, oscillations having a frequency approximately equal to that 
given by equation (23). Other frequencies, associated with wave travel in 
the cable, will of course also be present. 


Examples 


(1): Constant Dissipation Force; Maximum Initial Cable Tension 


Let us take as an example the following: 


Mg = 30,000 lb., mg = 3000 Ib., (m/M = 0.10), Vo, = 150 ft/sec., L = 50 ft., 

R = 150 ft., { = 250 ft., AE = 0.395 x 12 x 10 Ib., N=1. The frequency 

of cable tension oscillation given by equation (23) is f = 3.21 cycles/second. 
Assuming that the dissipation force is constant, we find from equation (6) 

that F, = 97,500 lb. It is assumed that the initial cable tension Ty is set at 


the maximum value consistent with the constant force, F,. Since N =1, 
To = Fe/2 = 48,750 Ib. 
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In order to carry out the solution we first determine x as a function of 
time using the simple energy relationship and the methods of the first section, 
and plot the result similar to that shown in Fig. 8a. Equation (22a) may then 
be solved for € in a step-wise manner by evaluating the constant 6 for each 
step and drawing the phase-trajectory. 

The initial conditions for the solution are 


If a time increment of 0.02 seconds is used the incremental angle in the 
phase-plane is given by 


p-At = 0.403 rad = 23° - 5', 


The phase-trajectory for € has been shown in Fig. 9, and the time- 
projection for €, in terms of cable tension T = AE €/) , in Fig. 10. The 
initial value of T is 48,750 lb. and the absolute maximum is 86,000 lb. The 
amplitude of the tension oscillation is approximately 48 per cent of the initial 
value. 

It should be noted in this example, and in the one that follows, that the 
mechanical advantage N of the pulley system at the arresting engine has been 
taken equal to one. This results in an abnormally large amplitude. It can be 
shown that if N equals 6 or 8, as may be the case in an arresting system 
on board ship, the amplitude is much smaller. 


(2): Constant Dissipation Force; Zero Initial Cable Tension 


The initial cable tension is zero, otherwise the constants are the same as 
in example (1). Since, in these examples, the energy absorber is assumed to 
be a dry-friction device no motion takes place in it until the force applied be- 
comes equal to the sliding friction force, that is when 2NT = Fe. (Actually, 
the cable tension must build up to a value somewhat greater than indicated 
due to the fact that the coefficient of static friction is greater thau the coef- 
ficient of sliding friction.) 

During the early stage of the arrestation the cable tension is given by the 
following, which has been obtained by setting s = 0 in equation (16): 


= Ke = # (24) 


Substituting expression (24) for T into equation (14) one obtains the following 
differential equation of motion for the aircraft 


oe 2 Lx 
x + x= Py 5 (25) 
+L 
where 
2AE 
gM 
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FIG. 9 PHASE-PLANE SOLUTION FOR € EXAMPLE (1) 
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Equation (25) may be solved readily by phase-plane methods by putting it in 


Lx 
the “6 form” where 8 = = = . The solution has been shown in 
242 
Fig. 11. From the resulting values of x and x it is then possible to deter- 
mine € and é from the following relations 


(26) 


and to determine the cable tension T from equation (24). 

At the instant when T reaches the sliding friction value, F,./2N, the energy 
absorber starts to slip and € and T are thereafter governed by equation (22a). 
The starting values of € ,€, and x to be used in the solution of equation (22a) 
are equal to the values of €, € and x determined from equations (24), (25), 
and (26) at the instant when T becomes equal to F;/2N. The phase-trajectory 
for x and x, after slipping starts, has also been shown in Fig. 11. 

Figures 12 and 13 show the phase-trajectory for € and € , and the time- 
projection of cable tension, for the initial stage of the motion before the 
energy absorber starts to slip as well as for the later stage during slipping. 

If the initial cable tension lies between zero and F,/2N the analysis fol- 
lows the same general procedure as in example (2), equation (25) governing 
until T becomes equal to F,./2N and equation (22a) governing thereafter. 


Experimental Investigation 


Experiments with a 1/10 length scale model, using a dry-friction strip 
brake as the energy absorber and having a mechanical advantage N = 2, have 
shown reasonably good agreement with the analysis.(4) Since the cable 
tension was recorded near the anchor end the effects of wave travel in the 
deck pendant were considerably attenuated by the reflections at the deck 
sheaves and by friction in the cable, and the recorded tension should there- 
fore agree well with the analysis. The results have been summarized below. 


Initial || Experimental Results Analytical Results| 
Initial || Xe 

Tension Number of| 
lin Cable 


A 
Tests prrsad | 


vergge 


1.28 to 1.31 | 1,30 
1.12 to 1.17} 1.14 


The frequency of oscillation in the cable tension, computed by equation (23) 
was 118 cycles/second, and determined from the experiments was 109 cycles/ 
second. 

It should not be concluded that the results of the model tests justify the 
neglect of the distributed nature of the cable mass. It is known from full 
scale, as well as model, tests that the maximum dynamic tensions in the 
deck pendant are affected to a considerable extent by wave travel. The 
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FIG. II PHASE-PLANE SOLUTION FOR x - EXAMPLE (2) 
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FIG. 13 CABLE TENSION — EXAMPLE (2) 
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effect of wave travel has been the subject of extensive analytical investigation 
by F. O. Ringleb, of full scale tests by the Naval Air Material Center and 
other agencies of the U. S. Government, and of model tests by the authors. 

A report on the model tests, including discussions of the scaling relation- 
ships, the choice of materials for the model cable, and the use of time- 
distorted models, will be available in the near future. 


Effect of Varying the Parameters - 
Constant Dissipation Function 


This section will indicate the effects of varying the initial cable tension, 
the mass concentration at the arresting engine, and the mechanical advantage 
in the cable and pulley system. It is assumed that the engine dissipation 
force F is constant and for N = 1 is equal to F, = 97,500 lb. The constants 
are Mg = 30,000 lb, V, = 150 ft./sec., L = 50 ft, J = 250 ft, R=3L, and 
AE = 0.395 x 12 x 108 Ib. 


a) Initial Cable Tension Varied; N = 1 


Fig. 14 shows the cable tension as a function of time for three values of T,, 
and for a mass ratio m/M = 0.10. We see that the maximum cable tension 
decreases as the initial cable tension is increased. 


b) Mass Concentration Varied; N = 1 


In Fig. 15 the cable tension is shown as a function of time for several 
values of the mass ratio and for the initial cable tension Tp = F,/2N = 
Fc/ 2 =T¢. The increase in maximum cable tension as the mass concentra- 
tion is increased is apparent in Fig. 15. Fig. 16 also shows the variation of 
maximum cable tension with initial cable tension and mass ratio. 


c) Mechanical Advantage Varied; Ty = Tc, m/M = 0.10 


In Fig. 17 the cable tension has been plotted as a function of time, for 
three values of the mechanical advantage. Fig. 18 shows that the ratio of 
maximum cable tension to constant cable tension decreases as the mechanical 
advantage is increased. It is also clear from Fig. 17 and equation (22b) that 


the frequency of oscillation of the cable tension varies directly as the mechan- 
ical advantage. 


Effect of Varying the Frontal Shape and 
the Rise-Time of the Dissipation Function 


A dissipation function that is constant from the start represents a condi- 
tion that does not normally occur in an arresting system and that is not 
necessarily desirable. However, it is fairly obvious from a simple energy 
analysis that the dissipation function F should be nearly constant over a 
large portion of the time of arrestation. Questions regarding the initial 
build-up or rise of the dissipation function do exist. Among the more im- 
portant are: Is the frontal shape of much consequence? Is the duration of 


the build-up or rise time of significance? The purpose of this section is to 
investigate these questions. 
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FIG. 1S CABLE TENSION MASS RATIO VARIED 
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The effect of varying the frontal-shape has been investigated for several 
different shape functions by calculating the response of the system to these 
functions. Fig. 19 shows the assumed functions and their mathematical ex- 
pressions. The constants of the system have already been given. The ratio 
of the rise-time, t; to the natural period T has been varied from 0 to 2.0 
(tT = 27/p, where p is given by eq. (22b) ). The mass ratio is m/M * 0.10 
and the mechanical advantage is N= 1. 

Fig. 20 shows the cable tension as a function of time for the versed-sine 
frontal shape and for various rise-times. The dotted curves indicate the 
static component, the maximum level being T, = Fc/2N = 48,750 lb. The 
cable tension curves show the expected oscillatory nature of the response. 
They also indicate considerable variation in the maximum cable tension as 
the rise-time varies. 

The variation is further illustrated in Fig. 21, where the spectra of maxi- 
mum cable tension versus rise-time are shown. It is evident that the effect 
of the rise-time is generally of much more importance than the effect of the 
frontal-shape function in determining the maximum cable tension. 


Comparison with a Simple Linear Oscillator 
and Step-Type Forcing Function 


The spectra in Fig. 21 are similar in general form to those for maximum 
displacement versus rise-time for an undamped linear oscillator starting 
from rest and acted upon by a step-type forcing function. Examples of the 
latter spectra and a schematic of the oscillator have been shown in Fig. 22.(9) 
The symbols are as follows: 


Ymax - Maximum displacement of mass M 


F 
Ys = 7 - static displacement of mass M 


F, - maximum constant value of forcing function 
k - spring constant of the oscillator 

ty - rise time of forcing function 

T = 2m /M/k - natural period of oscillation 


While these two systems have little direct relationship it is interesting to 
note that the relative importance of rise-time and of frontal-shape in deter- 
mining the maximum response is similar for both systems. In both types of 


systems for values of t;/7 greater than 1.5, the frontal shape is of little 
consequence. 


Effect of Varying the Parameters- 
Dissipation Function Not a Constant 


As before in the case of a constant dissipation function, we now investigate 
the effect of variation of the parameters when the frontal-shapes of the dissi- 
pation are not constant. Variations in both the mechanical advantage of the 
reeving system and the deck pendant length have been studied by calculating 
the response of the system using step-wise phase-plane methods. All of the 
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pertinent physical properties remain the same except where noted. 


~ a) Mechanical Advantage Varied 


The effect of varying the mechanical advantage for values of N ranging 
from 1 to 4 has been shown in Fig. 23 by the spectra of maximum cable ten- 
a sion versus rise-time. For this study the frontal shape has been limited to 
the versed-sine type with values of t,/7 ranging from 0 to 2. An increase 
in the mechanical advantage generally results in a decrease in maximum 
cable tension except in the vicinity of t;/T =1.0 where Tmax is nearly in- 


dependent of N. It may be noted for N = 4 that the optimum value of t;/rT 
is about 1.5. 


b) Deck Pendant Length Varied 


Fig. 24 shows the effect of varying the deck pendant length. The frontal- 
shape has been limited to the versed-sine case and the mechanical advantage 
is taken as N = 2. If the deck pendant length is reduced to zero the arresting 
system becomes a linear one for which an explicit solution is possible. This 
limiting case is discussed in the next section and is indicated by the dashed 
line in Fig. 24. A logical family of curves is formed by the spectra even 
though the curve for the limiting linear case is vastly different quantatively 
from those of the non-linear cases. 

The general trend is for an increase in the deck pendant length to result in 
a decrease in the maximum cable tension. However, the decrease is not very 
significant for values of L greater than 25 feet. From a comparison of the 
spectra shown in Figures 21, 23 and 24 it appears that the rise-time and the 
mechanical advantage are considerably more important than the deck pendant 
length in determining the maximum cable tension, except for very small 
(generally impractical) pendant lengths. 


* Analysis of the Limiting 
Linear Case (L = 0) 


The limiting linear case has no direct application to the usual type ar- 

. resting “gear.” However, the explicit nature of the linear solution makes it 
possible to explore qualitatively the effects of variation of the parameters in 
the non-linear cases. It also provides a limiting case, as already shown, for 
comparison with the cases discussed in the previous section in which the deck 

. pendant length is finite. Furthermore, the linear case is of interest in a more 
general study of systems for the arrestation of motion, not necessarily re- 
lated to aircraft arresting systems. 

The derivation of the differential equation for cable elongation and the ex- 
pressions for maximum cable tension, relating to the dissipation functions 
shown in Fig. 19, have been given in Appendix I. 


a) Effect of Varying the Frontal-Shape and Rise-Time 
of the Dissipation Function 


Fig. 25 shows the spectra of T,ax/T, versus ty/7 for various frontal 


shapes. The optimum value of ty/T for all shapes is about 0.5. All of the 


spectra tend toward the asymptotic value of 3.666 as indicated by the dashed 
line. 
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b) Effect of Varying the Mechanical Advantage 


The spectra of Tmax/Tc versus t,/7 for various mechanical advantages 
have been shown for the versed-sine front in Fig. 26. Again we see that the 
optimum value of t;/T is about 0.5. The asymptotic values of Tygx/Tc are 
3.666, 2.403, and 1.714, for N = 1, 2, and 4, respectively. The maximum 


cable tension shows a marked decrease as the mechanical advantage is in- 
creased. 


CONCLUSIONS 


It has been shown that the dynamical response of aircraft arresting sys- 
tems, treating them as “lumped parameter” systems and neglecting elastic 
wave travel in the cable, may be determined without great difficulty by step- 
wise linearization of the problem and by use of the phase-plane graphical 
method. This may be done for any form of the dissipation function (arresting 
engine force function). Furthermore, the analysis has been written so that 
cable elongation (or cable tension), which is the most important variable, is 
determined directly as the dependent variable rather than as the difference 
between two variables. 

The analytical expression (22c) for the step-wise constant 6 is a com- 
plicated function of aircraft displacement and of cable elongation. However, 
its evaluation presents no more than ordinary numerical difficulty since a 
good approximation for the aircraft displacement may be determined inde- 
pendently. 

The method is sufficiently rapid that a wide range of the parameters may 
be investigated, including the effects of varying the mechanical advantage, 
the deck pendant length, the masses, the initial cable tension, and the frontal 
shape and rise-time of the dissipation function. It should be pointed out that 
phase-plane methods may be used to advantage in the analysis of a wide 
variety of transient linear and non-linear problems in structural 
dynamics.(10,11,12) 

It is shown that there is a qualitative relationship between the response in 
cable tension of the non-linear arresting system and the response in displace 
ment of a simple linear oscillator acted upon by a “step-type” forcing func- 
tion. 

The following specific conclusions relate to the aircraft arresting system: 
The maximum of the “lumped parameter” oscillatory component of cable 
tension may be decreased by increasing the initial cable tension, by decreas- 
ing the moving mass associated with the arresting engine (energy absorber), 
by increasing the mechanical advantage in the cable reeving system at the me «4 
arresting engine, generally by increasing the length of the deck pendant, and “a q 
in general by increasing the length of time (rise-time) required for the dissi- 
pation force in the arresting engine to rise from the initial value to the maxi- -§ 
mum. The effect of varying the frontal-shape of the dissipation function is 4 
not of great consequence. i 


In the design of a particular arresting system, the deck pendant length and q 


the mechanical advantage in the cable reeving system are governed to a con- 3 
siderable extent by the space limitations. On board ship, for example, the ‘@ 
mechanical advantage may be 8 or 10, while on land it may be as low as 1 or zs 


2. Of the remaining parameters, those which appear to offer the greatest 
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possibility for reducing the oscillatory cable tension are an increase in the 


initial cable tension and a lengthening of the rise-time of the dissipation 
function. 
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APPENDIX I 


Derivation for the Limiting 
Linear Case (L = 0) 


If in Fig. 1 we let L equal zero the differential equations of motion may 
be written as follows: 


for aircraft: 
Mx = - 2T 
for the dissipation element: 


We can define 
and noting that 


x = NS+€ 
we obtain from equations (27) and (28) 


x= = Ns 


mé = - F +2NT (28) 4 

l 

t= = (x Ns) 

(29) 

(30) 

4 
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Then combining equations (29) and (30) we obtain 


, AE, m, x2 
=5 (31) 
Equation (31) could also have been obtained directly from equations (22) by 


letting L = 0. This differential equation has been solved for the various dissi- 
pation functions shown in Fig. 19. 


Expressions for Tmax, the maximum cable tension, occurring after the 
dissipation function has risen to a constant value, have been given below: 


Constant Slope Front 


T 
max _ 2 \ (32) 
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The asymptotic value of Tmax/Tc is given by: # 
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ABSTRACT 


Plastic Design Methods assume that local buckling of flanges and webs of 
WF-Beams will not occur during the formation of plastic hinges. Such severe 
conditions made the re-examination of the problem of plate buckling in the in- 
elastic range necessary. 

In contradiction to accepted opinions it was found that steel columns and 
plates can be compressed beyond the yield point and even into the strain 
hardening range without buckling. Theoretical results for the required geo- 
metric proportions are presented and comparison is made to experimental re- 
sults obtained from tests on model columns, angles and WF-beams. 

Furthermore, consideration is given to the problem of buckling between 
the elastic limit and the yield stress. 


INTRODUCTION 


2.1 Statement of Problem 


Present-day elastic design of steel structures is based on the concept of a 
specified safety factor against nominal yielding of the most stressed fibers. 
Concerning the detail design of flanges and webs of WF-beams and other plate 
elements, it is therefore sufficient that the yield stress can be reached with- 
out premature local buckling. However, newer methods of structural design, 
referred to as “Plastic Design”, “Collapse Design”, etc., propose to base the 
design on the actual load-carrying capacity of the structure. The working 


Note: Discussion open until September 1, 1958. To extend the closing date one month, 
a written request must be filed with the Executive Secretary, ASCE. Paper 1581 is 
part of the copyrighted Journal of the Engineering Mechanics Division, Proceedings 
of the American Society of Civil Engineers, Vol. 84, No. EM 2, April, 1958. 
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loads are determined as a specified percentage of the ultimate load. This 
ultimate load can only be realized if the members can undergo plastic defor- 
mations at a number of sections without local or lateral buckling, and a con- 
sequent fall-off in bending moment resistance. This process is generally re- 
ferred to as the formation of plastic hinges and redistribution of moments. 
It is therefore evident that Plastic Design imposes more severe conditions on 4 
the local buckling characteristics of the plate elements of structural members ££ : 
and requires a reappraisal of the problem of inelastic buckling. : 
Extensive literature* exists on the buckling of columns and plates com- 
pressed beyond the elastic limit of the material. However, it is generally ac- 
cepted that such elements made out of structural steel will invariably buckle 


once the yield stress is reached. The following statements from reference (1) 
may be quoted: 


On page 159: 


“To obtain a compressive load on a bar larger than the load produc- 
ing the yield-point stress, it is necessary to prevent the bar from 


lateral buckling at the yield point by applying some lateral con- 
straint.” 


And on page 385: 


“Experiments show that, when the compressive stress reaches the 


yield point of the material, .... the plate buckles for any value of 
the ratio b/h.” 


Essentially the same conclusions must be drawn from the equations for 
inelastic buckling presented in reference (2). Both the column and plate 
buckling equations—Eqs. (21) and (653)—contain the factor T = E,/E. When 
the yield stress is reached, the tangent modulus EF; reduces to zero and hence 
buckling seems to be unavoidable. 

At first glance, therefore, it may appear impossible to compress structur- 
al steel elements beyond the initial yield strain, € >, (see Fig. 1a) without 
buckling. This situation would invalidate one of the basic assumptions of 
plastic analysis and hence make it inapplicable to structural design. 

Fortunately this reasoning is in error. Whereas theories, as presented in 
references (1) and (2), give a fair account of buckling in the elastic range and 
in the inelastic range between the proportional limit and the initial yield point ‘ 
strain, € ¢, they are not applicable to cases beyond this range, A study of the 
stability of plates stressed into the strain-hardening range has been present- 


ed in reference (3). This paper constitutes a continuation and extension in- 
cluding pertinent test data. 


2.2 Inelastic Behavior of Structural Steel 


For a proper understanding of the new approach presented in this paper a 
short discussion of the inelastic behavior of structural steel is necessary. 
Fig. la shows the well-known stress-strain diagram of a coupon in tension or 
compression. In the elastic range the material exhibits homogeneous and 
isotropic behavior. Yielding commences rather abruptly without any signifi- 
cant indication of a proportion limit. At the yield stress, 0,, considerable 
straining takes place without an increase in stress such that the tangent 


*For a summary, references (1) and (2) may be consulted. 
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modulus, E;, seems to reduce to zero. However, it should be kept in mind 
that the strain, € , is an average strain, derived from the measurement of an 
elongation (or shortening) over a certain gage length. Actually there is no 
material within this gage length at a strain between € ¢ and €,. The mecha- 
nism of yielding is discontinuous, taking place in small slip bands by a sudden 
jump of strain from € gto €,. The slip bands form successively, starting at 

a “weak” point (inclusion, point of stress concentraction, etc.) and then spread 
out into the rest of the specimen.* 

Therefore, during yielding some of the material is still elastic while the 
yielded zones have reached the point of strain-hardening, €,. The material 
within the coupon is therefore heterogeneous. When all material has been 
strain-hardened, the stress starts to increase again. Once again, all of the 
material has identical physical properties and hence is homogeneous. How- 
ever, slip produced such changes that it is no longer isotropic—i.e., its 
properties are now direction dependent. If proper consideration is given to 
these physical facts the behavior of compression elements of structural steel 
in the yield and strain hardening range can be explained and predicted. 


Theoretical Investigation 


3.1 Inelastic Column Buckling 


Before treating the problem of plate buckling, it may be advantageous to 
first summarize the behavior of columns. The elastic buckling stress of a 
slender column is given by: 


(en (1) 
where 0. = elastic buckling stress 
E = modulus of elasticity 
< = effective column length 
r = radius of gyration 


In the inelastic range the theoretical stress at which bifurcation of equi- 
librium will occur (bending commences) is characterized by the tangent 
modulus stress‘2 


(2) 
Cn? 


where 


EF, = tangent modulus 


The maximum stress lies between Ot and the reduced modulus stress 


E, (3) 
Sr 


*For an informative description see reference (4) p. 297. 
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For a rectangular section the reduced modulus is given by(1) 


(\E (4) 


For design purposes the possible increase of the maximum stress beyond the 
tangent modulus stress is generally neglected. If the tangent modulus EF; is 
given, the behavior of a column can be predicted. 

On the basis of the simplified stress-strain diagram (Fig. 1b) buckling 
would be elastic up to the yield point stress, 04. The slenderness ratio (¢ /r) 
corresponding to initiation of yielding is obtained by inserting into Eq. (1) the 
values Oe = 0, and E. For smaller values of @/r buckling will not occur 
immediately once the yield strain €¢ is reached but yielding will commence, 
with the formation of slip bands somewhere along the length of the specimen. 
For a pin-ended column the worst situation will occur if yielding starts to 
spread from the center. The bending stiffness of the middle part is greatly 
reduced to E,I, Ep being the tangent modulus of the zone which has reached 
strain-hardening. The end-pieces are still elastic and hence their bending 
stiffness is given by EI (see Fig. 2, case a). Assuming that the column is 
continuously compressed such that buckling will occur without strain reversal, 
the buckling condition is given by the following transcendental equation:* 


E/E tan = | (5) 


where 
2 
P* = 
Ocy = buckling stress 


oO e = elastic buckling stress 


= parameter specifying extent of yielded zone as shown in Fig. 2, 
case a), 


For given values of E,,/E and ¢, values of ¥ can be derived by trial and 
error. As buckling occurs within the plastic range—i.e. the critical stress 
Ocr, is equal to the yield stress To the slenderness ratio ¢/r can be comput- 
ed: 


Ger = So = 960 = 
or: 
(t/r)* =¥ = (6) 
The shortening of the column, A 4, is given by: 


and the average strain at which buckling occurs: 


*See Ref. (1), p. 128, (Some changes in notation have been made). 
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Fig.2 Comparison of Theoretical Curves with 
Results of Model Column Tests 
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Ecr = at/et (1-25) +25€. 
Numerical results for Eq. (7) are shown in Fig. 2, curve (a), derived on the 


basis of the following material constants: 


E = 30,000 ksi 


E,= 800 ksi 
1.1 x 1073 
= 14x 1073 


" 


33 ksi 


The yield stress,o,, is reached for a slenderness ratio esr = 94.9, the 
corresponding strain being the yield strain € ¢- For smaller slenderness 
ratios, ¢/r, the buckling stress o;; remains constant at the yield stress 
level. However, buckling will only occur after the yielded zone at the center 
has spread over a length 22 ¢ resulting in a great reduction of the bending 
stiffness over this section. 

The average critical strain changes at first very little with decreasing 

€/r, but increases rapidly for values of @/r smaller than 30. Finally, for 
£/r = 15.4 the entire length is strain-hardened before buckling occurs. If 
@/r is still further reduced the critical stress will increase above the yield 
stress 0,. All material being strain-hardened, the tangent modulus E, for the 
entire length is given by the slope of the stress-strain curve in the strain- 
hardening range. 

Another extreme assumption would be that yielding spreads symmetrically 
from both ends. The corresponding curve (b) in Fig. 2 is considerably above 
the previously derived curve (a). 

The foregoing reasoning is substantiated by the results of 37 tests on small 
model columns which are also plotted in Fig. 2.* The average yield stress of 
the steel used was about 35 ksi. Unfortunately no records of the strain- 
hardening strain €,, and the strain-hardening modulus E,, are available. Al- 
most all test points lie between the two theoretically derived curves. The 

e actual yielding may occur anywhere between these two extreme assumptions. 

The experimental indications however are that strains larger than the yield 
strain € +, can be reached without buckling. For sufficiently small é/r values 
the point of strain-hardening, € 9, can be reached and even exceeded. Corre- 

° sponding test results will be given in a subsequent section. 

If theoretical predictions, based on Eqs. (1) and (2) and a tangent modulus 
Ey derived from a coupon stress strain diagram, are compared to the actual 
behavior of structural steel columns such as WF-sections, built-up members, 
etc., a considerable discrepancy is observed in the intermediate column range 
30 < &@/r <110, the theoretical predictions being always considerably too high. 
This situation generally has been attributed to initial imperfections, accidental 

end eccentricities, etc. Bleich introduced an “effective” coefficient, T = E;/E 
—Eq. (64) reference (2)—much below the actual T based on a coupon stress 


*The test results were obtained from Mr. P. C. Paris, instructor at Lehigh 
University. The tests were performed by him and Dr. T. A. Hunter at the 
University of Michigan under the direction of Prof. J. A. Van den Broek. 
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strain curve, in order to preserve the applicability of Eq. (2) to inelastic 
buckling. However recent theoretical and experimental investigations cleared 
up the situation very conclusively by demonstrating that this discrepancy was 
due to the presence of residual stresses.(5,6,7)* The proportional limit of an 
axially loaded column is reached when the sum of the applied stress P/A and 
the maximum compressive residual stress reaches yielding. The subsequent 
behavior of the column depends on the distribution of the residual stresses 
and the direction of buckling—i.e., weak or strong axis buckling. However for 
&€/r approaching about 15, the strain-hardening range is reached and the 
influence of residual stresses is completely wiped out. From there on the 
buckling load is governed by the tangent modulus, EF}, only. 


3.2 Inelastic Plate Buckling under Uniform Compression 


In the previous section the possibility of a column buckling at stresses 
above the yield stress,o,, was shown. It can therefore be assumed that the 
same may hold true for plates. The objective of the following analysis is the 
derivation of a plate buckling equation which is applicable to the strain-harden- 
ing range. The behavior of plates which buckle in the intermediate range be- 
tween the proportional limit (sum of applied and residual stress equal to yield 
stress) and the strain hardening range is governed by the presence and distri- 
bution of residual stresses. No direct solution of this problem has yet been 
developed. However, a reasonable transition curve will be proposed. 

During yielding the material changes its physical properties such that at 
strain-hardening the initially isotropic material has become orthotropic—i.e., 
the properties are direction dependent. In case of plane stress the stress and 


strain relations of such a material can be expressed in the following general 
form: 


normal strain 

shear strain 

normal stress 

shear stress 
= tangent modulus in x direction 
= tangent modulus in y direction 


= tangent shear modulus 


*Reference (6) contains a review of some of the experiments and earlier 
thinking that led the Column Research Council to take up the concept of Re- 


sidual Stresses and their influence on buckling and to stimulate research on 
this important topic. 
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Vs coefficient of dilatation for increase in Ox 


i. 2 coefficient of dilatation for increase in Oy 


If Eqs. (8) are valid for the entire plate cross-section, the expressions for 


the bending and twisting moments in terms of the transverse plate deflection 
w, become 


x 
= 


ox? dy? 
Jw ow 
My ($4 + (9) 
dw 
Mey 26, dX dy 
where 
t? 
t = thickness of plate 
E 
D, l-v 


The above expressions imply that buckling occurs without strain reversal. 
From these relations the following general expressions for the buckling 


strength of uniformly compressed, rectangular plates can be derived as shown 
in Reference (3): 


(1) Loaded edges x = 0 and x = @ hinged, unloaded edge y = 0 hinged 
and unloaded edge y = b free. (see Fig. 3) 


Ser = (+) (2 +G,} (10) 


where 
t = thickness of plate 
b ~~ = width of plate 
= 


length of plate or half-wave length 


For a long plate the first term can be neglected and the buckling stress be- 
comes: 
2 
(11) 


(2) Same as case 1, but unloaded edge y = 0 fixed. 


The minimum value of the buckling stress occurs when the length-to- 
width ratio is 


D 
|. 46 D, (12) 
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Fig.4 Plate Supported Along All Four Edges 
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Fig. 5 Dimensionless Representation of the Buckling 
Strength of Columns and Plates 
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Ger = {0.169 - 0.270 (Dy Dy.) + 1.712 Gel 


(13) 
Dyy Vy D, Dyx = Vy Dy 
(3) Loaded edges x = 0 and x =@ hinged, unloaded edges y = - d/2 


hinged. (See Fig. 4) The minimum value of the buckling stress 
is obtained for the following value of the length-to-width ratio: 


4) Dy 
(14) 
d | 


Gey = + Dy + Da 4G} as 


(4) Same as case 3, but unloaded edges y = ~ d/2 fixed. 
The minimum value of the buckling stress is obtained for 


Then 


= 0.66 (16) 
Then 
© (ty) { 4.554 (D,D, + 1.237(Dy+ Dy) + 


Several theories of plasticity are available for the determination of the 
moduli E,, Ey, vy, ¥, and Gy. The various theoretical predictions as applied 
to the plate buckling problem are summarized in Table 1. Obviously, there 
are significant discrepancies between these theories. From the point of view 
of the mathematical theory of plasticity the stress-strain relations should be 
of the incremental type, only the stress-strain law used by Handelman and 
Prager(12) satisfies this condition. All theories except the last one can, 

therefore, be discarded. 

Although the stress-strain law of Reference (12), with the second invariant 
of the deviatoric stress tensor as the loading function, has been verified ex- 
perimentally by combined compression and torsion tests on tubes of mild 
steel stressed into the strain-hardening range, (13) it cannot be a plied without 
modification to the plate buckling problem. Onat and Drucker‘ (4) pointed out 
that the Handelman-Prager theory is unable to predict the buckling of actual 
plates because of the influence of unavoidable initial imperfections. 

In Reference (3) it is shown how effective values of the moduli can be ob- 
tained from an incremental stress-strain relationship. From the stress- 
strain curve for the strain-hardening range of mild steel 


S-Ge 
= + k (18) 


en 
with 
E, = 900 ksi } 
¥ 
k= 21 
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m =2 


the following values of the moduli were found to be applicable:(3) 


4 
ve 


D, = 3,000 ksi 
Dy = 32,800 ksi 
Dyy = YyDy = ¥xDy = 8,100 ksi 


G, = 2,400 ksi 


The above developed expressions for the buckling strength of orthotropic 
plates can only be applied if all material is strained into the strain-hardening 
range (€ > €,). In the intermediate range transition curves can be developed 
as in the case of column buckling (see Fig. 2). 

Consider for instance a flange of length 24. The loaded edges x = 0 and 
x = 2@ are fixed, the unloaded edge y = 0 is hinged and the unloaded edge y = b 
is free.* For the simplified stress-strain curve of Fig. 1b the elastic so- 
lution is valid for 0.,.< 0, and obtained from Eq. (10) by substituting the 
elastic values of the moduli. Thus 


(LY (By +6 
b/ 12(I-*) \e (19) 

If the geometry of the specimen (given by b/t, b/g) is such that Ocy from Eq. 
(19) would exceed the yield stress, Oo, yielding will occur before the plate 
buckles, the stress remaining at the yield value ¢,,. In order to find a so- 
lution for this case it is necessary to assume a certain distribution of the 
strain-hardened zones, for example, yielding starts at the ends and spreads 
towards the middle. This assumption seems to be reasonable in view of the 
fixed-end restraints and has also been observed during tests. 

The middle section, being still elastic, is practically rigid compared with 
the yielded zones of length @. Assuming that only the latter will deform 


gives 


The corresponding average critical strain is 


Ecr = (1-5) + (21) 
Eqs. (20) and (21) are applicable to plate elements that are free of residual 
stresses (for example, annealed specimens). However, as delivered specimens 
contain in general residual stresses of considerable magnitude such that 
partial yielding will set in at an applied stress considerably below the yield 
stress. In the following a more realistic approach to the range 0, < d.,< % 
will be discussed. The elastic solution is only valid up to a limiting stress 
0,. The magnitude of g,, is determined such that the sum of the appliedstress 
op and the maximum residual compressive stress op equals the yield stress 


*This example has been chosen as a basis of comparison of the theory with 
the results of tests on plates with these boundary conditions presented in 
Chapter 4, 
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Oo: Hence the stress op corresponds to the “effective proportional limit” of 
the specimen. 


3.3 Buckling at Stresses Below the Yield Stress 


The elastic buckling stress, oo, of a perfectly plane plate of isotropic ma- 
terial, subjected to forces acting in its plane, is given by 


TE 
(6) 


where 
k = plate buckling coefficient 


This equation and the expression for the elastic buckling stress of a column, 
Eq. (1), can be written in dimensionless form as follows: 


G, | (23) 
— 
G. 
where for a column 
e Go 
(24) 
Tr IE 
and for a plate 
b | (26.(i-v 
a 
ait kE. (25) 


Eq. (23) is valid for values of o larger than a certain limiting value a@,, (See 
Fig. 5). The corresponding limiting stress is o,. As indicated before the 
sum of this limiting stress and the maximum residual stress equals the yield 
stress g,. From this point, (op/oo.%p), a transition curve must be followed 
to the point where the buckling stress equals the yield stress, (o¢;/a, = 1, 
Qo). If a specimen reaches this point all its material has been yielded and 


reached the strain-hardening range. The transition curve can generally be 
taken in the following form 


So S. ° (26) 

For a column the value of a, is relatively small and the influence ofstrain- 
hardening is, therefore, generally neglected by taking a, = 0. It has been 
shown by Huber, (5) that the shape of the transition curve for WF columns de- 
pends on the distribution of the residual stresses and the axis about which the 
column buckles (weak or strong axis). The magnitude of the residuals de- 
termines the value of o,/g,. Actual measurements of residual stresses in 
WF shapes have furnished rather high values. In order to be conservative, 
Op/%, for most WF shapes should be taken equal to 0.5. 

For WF Columns the influence of residual stresses is most pronounced if 
buckling takes place about the weak axis. The transition curve was found to 
be approximately a straight line (n = 1, @,, = 0). When the column fails by 
bending about the strong axis a quadratic parabola, (n = 2, a, = 0) is a good 
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approximation of the transition curve. 

For plate buckling the influence of residual stresses will be less severe. 
Indeed, the buckling stress is not only dependent on the modulus Ey, in the di- 
rection of the compressive loads, but also on the modulus E,, in the transverse 
direction and the shearing modulus G; (see Eqs. (8), (10 to 17) and (18)). 
Yielding affects the latter to a much lesser degree than Ey. This influence 
can be taken into account through a higher value of n. However, the transition 
curve given by Eq. (26) can at the most be tangent to the elastic solution given 
by Eq. (23). From this condition the limiting value of n is found: 


2(X.- %e) 
max (<? \) 


It is suggested that this value be used for all cases of plate buckling. 

Finally, the problem remains to determine the values of ap andao. Eq. 
(23) gives the value of ay by substituting 0, = 0p, where dp is the effective 
proportional limit. For structural wide-flange shapes it is conservative to 
take = 0.5. Thus a@ 2. Furthermore, the following values of a, 
are determined from Eqs. (2), (11), (13), (15) and (17) respectively by substi- 
tuting the given values of the moduli. In parenthesis the corresponding values 
of ¢/r, b/t or d/t are given for 0, = 36 ksi. 


(27) 


Columns = 0.173 = 15.7) 
Long hinged flanges: ot, = 0.455 (b/t = 8.15) 
Fixed flanges : Bo = 0.461 (b/t = 14.3) 
Hinged webs > @o = 0.588 (d/t = 32.3) 
Fixed webs > @o = 0.579 (d/t = 42.0) 


The above table suggests that the values of a, depend only on the type of 
compression element (column, flange or web plate) and are nearly independent 
of the amount of restraint. The problem can, therefore, be somewhat simpli- 
fied by taking for columns a, = 0.17, for flanges (plates with one free edge) 
G@, = 0.46, and for webs (plates supported along all four edges) a, = 0.58. 


These values are used subsequently in Fig. 23 where comparison to test re- 
sults is made. 


3.4 Inelastic Plate Buckling under Combined Bending and Compression 


An extension of the previous considerations to cases of plates subjected to 
bending and axial load is of practical importance. The web of a WF-beam 
subjected to a bending moment M and axial load P presents such a case. De- 
pending on the ratio of P/Py, Py being the yield load of the axially loaded 
member and equal to o,A, the neutral axis may lie inside or outside of the 
web. 

First of all, the plate buckling coefficient k must be determined. The mini- 
mum values of k for a stress distribution of a fully plastified section are 
shown in Fig. 6. In this figure, d corresponds to the depth of the web and y,, 
fixes the position of the neutral axis. The values were determined by equating 
the work of the external forces to the dissipation of the internal energy at the 
moment of buckling, using the moduli given in Chapter 3.2. 

Secondly, the strain €.., at which buckling of the web plate occurs and the 
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corresponding depth to thickness ratio d/t is of importance. For an expedient 
solution of this rather involved problem recourse to subsequently described 
test results is made. Fig. 7 shows experimentally determined values of a as 
a function of the critical strain €¢y for uniformly compressed webs. The 
parameter a is defined by Eq. (25)—with depth d instead of width b: 


12. (1-v) (28) 
wt | kE 


and is hence proportional to the d/t ratio. Turning now to the problem at 
hand, the maximum strain of the compression flange is taken as €,,. As the 
strain varies linearly over the depth of the section the mean strain of the 


compression zone in the web is € ,,/2. Using this latter value in entering Fig. 


7 the following values of @ are found: 
€m/€¢ = 12: a@ = 0.58 
€m/e¢= 8: = 0.60 
4: & = 0.69 
The resulting curves for the critical d/t ratios are plotted in Fig. 8 for the 


specific values indicated. 
Values of d/t corresponding to 0, different from go, = 33 ksi are obtained 


by multiplying d/t byy 33 ksi/ao. 


Experimental Investigation 


4.1 Compression Tests on Short Rectangular-Columns 


Contrary to commonly accepted opinion, it was stated above that buckling 
of columns at stresses above the yield stress should be possible without 
special precautions to prevent buckling at the yield stress. In order to verify 
this statement experimentally, a number of compression tests on short model 
columns of rectangular cross-section were carried out. The dimensions of 
the specimens, cut from the flange of an 8WF40 section, (ASTM A-7 steel) 
were chosen such that buckling in the strain-hardening range could be expect- 
ed (see Table 2). The test set up is shown in Fig. 9. The columns were 
placed flat-ended in a hydraulic testing machine. During the tests simultane- 
ous load, strain and lateral deflection readings were taken. 

The maximum loads divided by the area are plotted in Fig. 11, and com- 
pared with the tangent-modulus and reduced modulus loads. Confirmation is 
obtained of the prediction that the maximum load lies between the tangent and 
reduced-modulus load. The behavior of the columns is illustrated in Fig. 12 
showing the lateral deflection of the column as a function of the strain. It is 
obvious from these curves that the columns can reach the strain-hardening 
range without deflecting excessively. The critical strains corresponding to 
the tangent-modulus loads are indicated by arrows. Indeed, at these points 
the deflections start to increase at an accelerated rate, confirming once again 
that the tangent modulus load is the load at which bending commences. 


4.2 Compression Tests on Angles 


The next phase in the experimental program was the verification of the 
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Fig. G Plate Coefficient of Webs of Fully Plastified 
Wide-Flange Sections Subjected to Axial 
Load and Moment. 
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Fig. 7 Experimentally Determined Values of a 
as a Function of the Critical Strain of 
Uniformly Compressed Webs. 
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TABLE 2 


Dimensions of Short Columns 


b = width 
t = thickness 
L = length of specimen 


C/r = = effective slenderness ratio (test 


conditions simulate fixed ends) 


r "as" radius of gyration 
1 


we 
Be 
Specimen] in. in. in. 
+ 
C2 0.745, 0.544] 2.00 4 
Ch 0.745) 0.543] 3.65 11.6 
cS 745) 0.545] 4.33 | 13.9 
C6 0.745) 0.545] 4.80 15.3 
7 
C7 0.750] 0.527] 5.65 18.9 
0.750[0.527] 6.90 | 23.4 
7 
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Fig. 9. Short Column Test 
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Fig. 10 Angle Test 
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TABLE 3 


TABLE 


Results of Angle Coupon Tests 


Strain at 
Strain- 
hardening 


Type of 


Loading 


tension 


" 


As- 


elivered 


All coupons tested in Baldwin 60,000 1b. Hydraulic 


Machine. Valve opening corresponding to tes 


speed of 1 micro-in/in per second in the elastic 


range. 


7 
Dimensions of Angle Specimens ae 
ength idth | Thickness 
Specimen 2e b b/t | 2e/t Material 
fs 
in) in) (in) 
A-21 25.0 | 4.87 0.383 12.70] 5.14 ASTM 
A-22 25.0 | 4.79 0. 361 12.60} 5.21 A-7 ae 
A-31 i7:9 0.370 8.85] 5.48 Année 
A-32 17.9 3.26 0.374 5.79] 5.46 ealed- 
| 
12.5 | 2.31 0.377 6.13] 5.41 | 
A-42 12.5 26 34 0.371 6.36 De 35 
@ A-33 17.5 | 3.30 0.378 8.73] 5.30 ASTM i 
nee A-51 21.2 4.07 0. 380 10.70] 5.21 So 
As- 
fie delivered 
4 
Coupon] Material | Stress hardening 
E,x10° Eo ksi 
2 
Cl 37.6 14.3 906 compression 
2} " Be 
C2 Ann- 35.1 14.6 
Tl ealed 36.8 18.5 903 
3 
QLo 
T2 38.2 20.0 869 " or 
: T3 | 40.2 9.5 577 
Th 39.9 15.6 1,160 
4 
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Fig. 13 Results of Angle Compression Tests 
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inelastic plate buckling theory. For this purpose a number of compression 
tests on angles were carried out. When buckling torsionally, the flanges of 
an equilateral angle act as two plates, hinged along the heel and free along the 
outstanding edge. The loaded ends of the angle columns were fixed in the 
testing machine against rotation. Fig. 10 shows one of the specimens after 
testing. The dimensions of all specimens are given in Table 3. The speci- 
mens were machined from a 6 x 6 x 3/8 (14.9 lb) angle, ASTM A-7 steel, by 
reducing the width of the leg to the desired dimension b. Table 3 indicates 
that six specimens were tested after complete annealing and two specimens 
were in the “as-delivered” state, Material properties as obtained from 
coupon tests are summarized in Table 4. 

Fig. 13 shows the obtained stress vs. average axial strain curves. Further- 
more, the rotations of the center section were measured. They are plotted 
as a function of the average strain in Fig. 14. In both figures the experi- 
mentally determined critical strains are indicated by arrows. The critical 
strain is defined as the strain at which the rotation starts to increase more 
rapidly than it did initially. Specimens A41 and A42 did not buckle torsionally 
but failed by bending about the weak axis. Buckling occurred in the strain- 
hardening range, thus proving once more that buckling does not necessarily 
take place at the yield stress. The results of all tests on angles are summa- 
rized in Table 5. 

The theoretical solution for the strain-hardening range is given by Eq. (10) 
and for the yielding range by Eq. (20) and Eq. (21). Fig. 15a shows a plot of 
¢ (parameter defining extent of yielded zones) vs. b/t as obtained from Eq. 
(20). As elastic deformations have been neglected in the derivation of this 
equation b/t =-o for £ = 0 (rigid-plastic solution). From Eq. (19), however, 
it follows that for ¢ = 0, b/t = 20.7. Knowing the rigid-plastic solution and the 
point for £ = 0, corresponding to the limit of the elastic solutions, the elastic- 
plastic solution is sketched as a dotted line in Fig. 15a. The critical strain 
can then be determined from Eq. (21). 

The obtained € ., vs. b/t curves are compared with test results in Fig. 
15b. The theoretical curve describes adequately the behavior of the test 
specimens. 


4.3 Compression and Bending Tests on Wide- Flange Shapes 
Six wide-flange shapes were tested under two loading conditions: 


(a) Axial compression (Tests D1, D2, D3, D4, D5, D6) 
(b) Pure bending (Tests B1, B2, B3, B4, B5, 6) 


The dimensions of all specimens are given in Table 6. Material properties 
obtained from coupon tests are summarized in Table 7. The compression 
specimens were placed flat-ended in a testing machine, the length being divid- 
ed in three gage lengths over which the change in length was measured direct- 
ly with 0.0001" Ames dials. The set-up of the compression tests is shown in 
Fig. 21a. The bending specimens were loaded symmetrically by two concen- 
trated loads (see Fig. 21b) and the dimensions were chosen such that shear 
failure was avoided. Again changes in length of the flanges were measured. 
The resulting P/A vs. € ay curves for the compressions tests and M/Z vs. 
€ ay curves for the bending tests are plotted in Figs. 16 and 17 respectively, 
where 


P =compressive load 
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= area of cross-section 
M = moment 
Z =plastic section modulus 


€ay = average strain at center of compressed flange. 


= Furthermore, with the aid of a dial gage lateral deflection measurements 
were taken along the edges of the flanges and the center of the webs (see Fig. 
18). For the bending tests the lateral rotation was measured at mid-span of 
the beam. The loading points were supported against lateral rotation. The 
observations of web deflection, flange buckling and lateral rotation from test 
B3 are plotted in Fig. 18 as a typical example. From these curves the critical 
strains are obtained as indicated by arrows in Figs. 16 and 17. As mentioned 
above, the critical strain is defined as the average strain at which the de- 
flection starts to increase more rapidly than it did initially. The results of 
hy all tests are summarized in Table 8. 

as The results of those tests where flange buckling was predominant are 
shown in Fig. 19. The results of tests D4 and D6 are omitted because web 
buckling occurred first and caused premature flange buckling. Specimen B4 
did not develop a major flange buckle but failed by lateral-torsional buckling. 
This result is therefore eliminated from Fig. 19. The critical strains are 
plotted as a function of b/t and are compared with theoretical solutions* for 
different values of the coefficient of restraint,8. This coefficient represents 
the restraints provided by the web to the flange. In case of elastic restraint, 
where ~W = moment per unit length required for a unit rotation, the coefficient 
& of restraint becomes 

" From the tests results it is seen that 8 has a value of the order of 0.01. 

= For the cases where web buckling occurred first (Tests D4 and D6) or 

a simultaneous with flange buckling (Test D2) the results are plotted in Fig. 20 
ig and are compared with the theoretical solution given by Eqs. (15) and (17). 

ry It is concluded that the theoretical curves give a good description of the 
actual behavior of flanges and webs. 

Some typical examples of wide-flange specimens after testing are shown in 
Fig. 22. Specimen D6 failed by web buckling after just reaching the yield point 
stress. The yielded material is concentrated at the top of the specimen. 
Simultaneous web and flange buckling caused specimen D2 to fail after just 
being completely strain-hardened while specimen D5 failed by flange buckling 
in the strain-hardening range. Of the bending specimens B2 and B6 failed be- 


Ai cause of flange buckling but B4 buckled laterally without developing a flange 
buckle, 


4.4. Summary of Test Results 


All test results are summarized in Fig. 23 showing 0,,/¢, as a function of 
the dimensionless parameter a. The figure clearly illustrates the difference 
in the behavior of columns, flanges and webs. The values of a for which the 
strain-hardening range can be reached are confirmed by the test results and 


*The theoretical curves are taken from Reference 3. 
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Fig. 2la Wide-Flange Compression Test (D3) 


Fig. 21b Wide Flange Bending Test (B1) 
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Fig. 22 Typical Examples of Wide- Flange Specimens After Testing 
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can be taken as follows: 

Column Buckling a, = 0.17 

Flange Buckling O, = 0.46 

Web Buckling Gt, = 0.58 


For flanges and webs, the shape of the transition curves (the plotted curves 
correspond to the maximum value of the exponent n, Eq. 27)) should be con- 
sidered as reasonable interpolations in the absence of sufficient test results 
in this particular region. Concerning the transition curve for the columns 
considerable experimental and theoretical data are available. (5) 


Recommendations for the Geometry of Wide- Flange Shapes in Plastic Design 


5.1 Flanges 


The following recommendations are directly based on the above presented 
theoretical and experimental investigation. Flanges of WF shapes can be 
strained into the strain-hardening range if a < 0.46. In general this is a 
sufficient requirement for the development of plastic hinges. Furthermore, 
buckling is then not accompanied by a sudden drop of the ioad (see results of 
tests A31-32-33, D2 and B2 in Figs. 13, 16 and 17 respectively). Consequent- 
ly the following recommendation can be given: 


The ratio of the outstanding width to the thickness of flanges, b/t, 
shall not exceed the following values depending on the value of the 
yield stress, 


for 0, = 33 ksi: b/t < 8.7 and 
for 0, = 36 ksi: b/t < 8.3. 


A 


5.2 Webs 
(a) Uniformly Compressed Webs: 


The results of compression tests on WF shapes failing by web buckling 
showed that for o& = 0.77 (Test D6) the yield stress can just be reached and 
for a = 0.57 (Test D2) the section can be uniformly compressed up to the 
strain-hardening range. For a section subjected to an axial load only, the 
first attainment of the yield stress is generally a sufficient requirement. In 
this case, therefore, the following recommendation can be given: 


If a section is subjected to an axial load only, the ratio of the 
distance between the center planes of the flanges over the thick- 
ness of the web, d/t, shall not exceed the following values de- 
pending on the magnitude of the yield stress o,: 


for 09 = 33: d/t < 44, and 
for 09 = 36: d/t < 42. 
(b) Webs in Compression and Bending: 


WF-Members subjected to axial. load P and bending moment M may be 
required to deform plastically. Obviously the critical depth/thickness ratio, 
d/t, depends not only on the stress distribution but also on the required 
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maximum strain €,, of the compression flange. Curves relating the critical 
d/t to the stress distribution (parameter P/g,Ay) and a given maximum 
strain (parameter € p)/€ ¢) were derived in Chapter 3.4 and are shown in Fig. 
8. Based on these results the following recommendation is made: 


If a WF section is subjected to axial load and bending moment the 
d/t ratio shall not exceed the values given in Fig. 8. Except in 


cases requiring large rotations the curve for € m/e f¢ = 4 can be 
used. 


Values of d/t for a yield stress 0, different from 33 ksi may be obtained by 


multiplying the ratio by 733 ksi/o,. From the curves of Fig. 8 three specific 
examples of P/P,, vs. allowable d/t ratios were investigated and are shown 

in Fig. 24. It can be seen that the 12JR11.8 section, for example, is not 
recommended for use in plastic design, the 12B14 section can be used if 
P/P,,< 0.08 and the 36WF108 if P/P, < 0.135. 


Summary 


Contrarily to commonly accepted opinions columns and plates of structural 
steel can be compressed beyond the yield stress level and even into the strain- 
hardening range, provided certain geometric conditions are met. After short- 
ly discussing the discontinuous yielding process of structural steel the values 
of these geometric conditions as slenderness ratio @/r for columns and width- 
thickness ratio for plates are derived theoretically. The interesting obser- 
vation is made that plates reach the strain-hardening range with a relatively 
smaller reduction in the d/t-ratio than the corresponding reduction in @/r for 
columns. 

The findings are substantiated by a series of tests on rectangular bars in 
axial compression, on angles and WF sections in compression and WF 
members in bending. 

The results are used to specify the geometric cross-section proportions of 
WF beams such that these members may develop large plastic deformations 
without local buckling and a consequent fall-off in load. These later re- 


quirements are essential for a successful application of Plastic Design 
Methods. 
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Nmax = maximum value of n defined by Eq. (27) 
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width of plate supported along all four edges 


modulus of elasticity 


strain-hardening modulus 


reduced modulus 


secant modulus 
= tangent modulus 


= tangent modulus in x-direction 


tangent modulus in y-direction 


modulus of elasticity in shear 
= tangent modulus in shear 
moment of inertia 


coefficient in Eq. (18) 


plate buckling coefficient 

= length of column and WF compression specimen 

= length of WF specimen subjected to pure bending 

= effective length of column 

= half-wave length of buckled shape 

= bending moment 

= bending moment per unit width of plate in x-direction 
= bending moment per unit width of plate in y-direction 
= exponent in Eq. (18) 

= exponent in Eq. (26) 


= axial load 

= 0,A = yield load 
= radius of gyration 
= plate thickness 

= coordinate 


= coordinate 
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= distance between neutral axis and compression flange 
= plate deflection 
= dimensionless coefficient defined by Eq. (24) and Eq. (25) 
value of @ at point of strain-hardening 
coefficient of restraint 
shearing strain 
normal strain 
average strain 
= critical strain 
yield strain 
maximum strain of compression flange 
strain at beginning of strain-hardening 
coefficient determining extent of yielding 
angle of lateral rotation 
Poisson’s ratio 
coefficient of dilatation for stress increment in x-direction 
coefficient of dilatation for stress increment in y-direction 
normal stress 
critical (buckling) stress 
elastic buckling stress 
stress corresponding to limit of validity of elastic solution 
yield stress 
residual stress 
reduced-modulus stress 
tangent-modulus stress 
shear stress 
E;/E 


edge moment per unit length to produce unit rotation of edge 
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EFFECT OF DEFLECTION ON LATERAL BUCKLING STRENGTH@ 
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SYNOPSIS 


The effect on lateral buckling strength of the deflection or curvature that 
occurs in an initially straight beam or beam-column prior to buckling is in- 
vestigated. The results of an analytical study of beams and an experimental 
investigation of beam-columns are presented. The effect of deflection on the 
critical moment is found to be larger when the member has ends simply sup- 
ported in the lateral direction than when lateral restraint is provided. It is 
concluded that in the preparation of design rules for beams the effect of de- 
flection need not be considered and that, in the case of beam columns, the only 


effect of deflection that need be taken into account is the increase in bending 
moment caused by the deflection. 


INTRODUCTION 


In the derivation of expressions for the critical moment for lateral buckling 
of beams, the assumption is usually made that the deflections in the plane of 
the applied loads are small enough that they may be neglected. Some investi- 
gators have made the equivalent assumption that the beam has an initial 
camber just sufficient to nullify the deflection at the buckling load. 

In 1899, A. G. M. Michell (1) presented an expression for the lateral buckling 
strength of beams that included part of the effect of the deflection of the beam 
axis just prior to buckling. An assumption that limits the generality of 
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Michell’s solution is that the torsional rigidity is negligibly small compared to 
the principal flexural rigidity. In 1952, Pattersson, 2) obtained a more gener- 
al solution for the critical moment of a rectangular beam that reduces to 
Michell’s solution when the necessary assumption is made regarding the 
torsional rigidity. 

The effect of the deflection of an initially straight beam was taken into ac- 
count in a study of proposed British specifications for the design of beams. (3) 
For some of the beam cross sections considered in that investigation, this 
effect increased the calculated buckling strength by as much as 15 per cent. 

Numerous tests of the lateral buckling strength of beams and beam-columns 
have shown excellent agreement between theoretical buckling strengths and 
test results.(4,6) It has been pointed out that in obtaining this close 
agreement, (5 no account was taken of the effect of the deflection as analyzed 
by Michell. However, the beams and beam-columns in question had ends fixed 
against rotation in the lateral direction, whereas Michell’s solution, as well 
as the other theoretical treatments of this problem, applied only to beams with 
laterally hinged ends. The effect of the deflection is appreciably different in 
the two cases. To the writers’ knowledge there has been no previous study of 
the effect of deflection or curvature on buckling strength of beams with later- 
ally restrained ends. Therefore, a brief study of this problem was made to 


determine whether or not this factor should be taken into account in design 
rules. 


Discussion 


The deflection that occurs before lateral buckling in an initially straight 
beam subjected to uniform bending, causes the axis of the beam to become in- 
clined to its original position and, hence, changes the summation of forces on 
cross sections normal to the axis. This introduces modifications in the differ- 
ential equations of equilibrium, that account for the increased critical moment 
resulting from deflection found in theoretical investigations of beams with 
laterally simply supported ends.(1,2) However, one of the changes in the 
differential equations of equilibrium resulting from deflection appears only if 
the ends of the beam are restrained; the inclination of the cross sections which 
accompanies the deflection introduces a component of the lateral restraining 
moment that tends to twist the beam and thus weaken its resistance to lateral 
buckling. 

When the member under consideration is a column loaded eccentrically in 
the plane of its web, the critical moment for lateral buckling is affected in the 
same way as previously described for a beam. It is this effect of the de- 
flection on the critical moment that is the subject of this investigation. In ad- 
dition, the deflection increases the moment arm of the end load, which causes 
an increase in bending moment. This second effect of deflection has been 
taken into account in other investigations of eccentrically loaded columns. (4,6) 

This investigation was divided into two parts. The first part was an ana- 
lytical study of the buckling strength of I-beams under uniform bending 
moment with ends hinged in the lateral direction, and of deep rectangular 
beams with the ends fixed in the lateral direction. This part of the investi- 
gation was not extended to I-beams with ends fixed against lateral rotation. 
However, the study of rectangular beams serves to give some indication of 


the effect of deflection and the end conditions on the theoretical buckling 
strength of I-beams. 
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The second part of the investigation involved tests on eccentrically loaded 
I-section columns with ends simply supported about both the principal axes. 
These tests were compared with the results of tests on similar members hav- 
ing ends fixed in the lateral direction.(6) This comparison helped to indicate 
whether the difference in the effects of deflection caused by varying conditions 


of lateral end fixity applies to this type of member as well as to the rectangu- 
lar beams for which the analytical study was made. 


A. Theoretical Analysis 


If the effect of deflection on the equilibrium equations is taken into account, 
the critical moment for a narrow, rectangular beam with ends simply support- 
ed in the lateral direction and subjected to a uniform bending moment in the 
plane of its greatest stiffness is given by the following equation: (2) 


where M' = critical bending moment in the plane of greatest stiffness of the 
member under pure bending, in.-lb, 


L = laterally unsupported length of the member, in., 
B, = El; = least flexural rigidity of the member, lb-in.2, 
Bo = Elg = greatest flexural rigidity of the member, lb-in.2, 
C = JG = torsional rigidity of the member, Ib-in.2, 
E = modulus of elasticity, psi, 


I, = moment of inertia of the member about the axis in the plane of 
the applied bending moment, in.4, 


Ig = moment of inertia of the member about the axis perpendicular to 
the plane of the applied bending moment, in.4, 
J = torsion constant, in.4 


G = modulus of rigidity, psi. 


The quantity (1 - C/Bg) in the above equation does not appear in the original 
solution for the case of the narrow, rectangular beam given by Michell. (1) 
This is because one of the terms that appears in the equilibrium equations was 
omitted in Michell’s derivation, as noted in the Appendix. Consideration of 
the effects of deflection caused the presence in Eq. (1) of the factor: 


V (- gt) (- &) 


In the Appendix the following equation is derived for the critical bending 
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moment of a narrow, rectangular beam under uniform bending moment with 
the ends fixed against lateral rotation: 


(- 


By comparing this result with Eq. (1), it is seen that the correction factor has 
been modified by an additional factor of (1 - B;/Bg). 

In Fig. 1 are plotted the correction factors by which the critical moment of 
the deep rectangular beam is multiplied to account for the effects of deflection 
for the cases of laterally hinged and laterally fixed ends for various depth/ 
width ratios. It will be noted that the factor for laterally hinged ends is sig- 
nificantly greater than unity (10 per cent or more) for the smaller ratios of 
depth/width but approaches 1.0 rather rapidly for deeper beams. The factor 
for fixed ends is not significantly different from 1.0 for any depth/width ratio. 

It is also shown in the Appendix that for an I-beam loaded with equal end 
moments in the plane of its web and with ends simply supported laterally, the 
critical moment is increased by the same factor that applies to rectangular 
beams if deflection is considered. Thus, the effect of deflection on a laterally 
hinged member is the same, whether the member has a deep rectangular 
cross-section or an I-section. The theoretical analysis of the effects of de- 
flection on the critical moment of I-beams with ends laterally fixed has not 
been completed. For sufficiently slender I-beams, however, the effects of 
restraint of warping can be neglected, and Eq. (2) would apply. In this case, 
the correction factor for deflection would be small, as is the case for rectangu- 
lar beams with laterally restrained ends. 


M! 


B. Experimental Investigation 


A series of tests on eccentrically loaded I-section columns was conducted 
to determine whether the effects of deflection indicated by the theoretical in- 
vestigation apply to these members. Since a number of tests of eccentrically 
loaded I-section columns had previously been made with the column ends fixed . 
laterally,\4,6) the tests in this investigation were made only with ends simply 
supported laterally. 

The specimens were 55-in. long extruded H-section members of 2014-T6 
aluminum alloy. The cross section dimensions, section properties, and me- -§ 
chanical properties are given in Fig. 2. The ends of the column specimens 
were first finished flat and parallel by turning on an arbor in a lathe. The 
corners of the flanges were then machined off at 45-degree angles in a shaper, 
leaving a 1/32-in. wide flat surface at the center of the web. Thisflat surface 
was coated with Prussian Blue and the corners of the flat surface rounded off 
with a file, leaving a line of Prussian Blue in the center to insure that straight- 
ness and parallelism of the knife edges were maintained. This procedure re- 
sulted in a 90-degree knife edge, with about a 1/32-in. radius, at each end of 
the specimen as can be seen in Fig. 3. These knife edges were supported in 
120-degree V-grooves machined in steel bearing blocks to permit lateral ro- 
tation at the ends. Rotation in the plane parallel to the web was permitted by 
mounting the bearing blocks on trunnion-type ball bearing heads of 10,000-lb 
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capacity. Since the center of rotation of the trunnion heads is in the plane of 
their platens, the length between supports in the plane of the web was greater 
than the specimen length by twice the thickness of the steel bearing blocks. 
This difference in length was 1-1/2 inches, giving an overall length of 56-1/2 
in. in the plane of the web compared to a length of 55 inches in the plane 
perpendicular to the web. 

Electrical resistance strain gages of the SR-4 type were used on the inside 
face at each flange corner at mid-height of the column specimens. Lateral 
deflection readings were taken with a wire and mirror scale arrangement. 

As indicated in Table I, tests were made on four specimens at nominal 
eccentricities of 0, 1/8 in., 1/2 in., 3/4 in. and 7/8 in. The specimens with 0, 
1/8-in., 1/2-in. and 3/4-in. eccentricities were tested with the knife edges 
dry, and specimens with 0, 3/4-in. and 7/8-in. eccentricities were tested with 
knife edges lubricated with a suspension of graphite in grease. Those speci- 
mens which were tested at more than one eccentricity have letter suffixes in 
the specimen numbers, the letters indicating the order of testing. All loads 
were applied in the plane of the web of the specimens. 

The specimens were loaded in a 40,000-lb capacity Amsler testing machine, 
using the zero to 10,000-lb range. The critical axial load was determined by 
tests on the two specimens with the zero nominal eccentricity. The critical 


load can be expressed by the Euler formula for buckling in the elastic range 
as: 


TPE 
or 


(KeL2)* 


whichever is smaller, where the subscripts 1 and 2 refer to buckling in the 
directions perpendicular and parallel to the web, respectively, and K is the 
effective length coefficient, depending on end conditions. Both the specimens 
tested concentrically failed by lateral buckling. The values of K; correspond- 
ing to the measured critical loads were 0.989 for the dry knife edges and 
0.997 for the lubricated knife edges. On the bases of previous tests with the 
trunnion heads, (4) it was assumed that the effective length factor Kg for buck- 
ling in the plane of the web was 0.99, giving an effective length of .99 x 56.5 
in. Using this information, the quantity P" was calculated to be 27,980 lb. 

The specimens with eccentrically applied loads were then tested. Knowing 
P' and P", it was possible to calculate the eccentricities from the strain data 
by means of the secant formula: 


A= e (sec VX -1) (3) 


where A is the deflection of the column, inches 
e is the eccentricity of the applied load, inches 
a@ is the ratio P/P". 


This expression can be put in the form: 


(M-Mo)c ec T 17 
P sec > Po see V% | 
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where M = Pe = bending moment, in.-lb 
Mo = Poe = initial bending moment, in.-lb 
P = end load, lb 
P,, = initial end load for strain measurements, lb 
c = distance of strain gages from neutral axis, in. 
Qo = Po /P" 


The left side of Eq. (4) represents the measured bending strains. The 
value of e can be determined graphically from this equation as in Fig. 4. The 
values of e thus determined were used to obtain the value of M at failure. The 
specimen with zero nominal eccentricity and dry knife edges showed no 
measurable eccentricity, but the specimen tested nominally concentrically 
with lubricated knife edges had an actual eccentricity of 0.008 in. This was 
not sufficient to affect the experimental value of the critical load, since the 
calculated P/P' ratio at failure for this eccentricity is 0.99999. 

The strength of beam-columns is often computed by means of interaction 
equations that relate the strength under combined load to the critical loads 
under axial load or end moments acting separately. Such interaction curves 
are shown with the test data in Fig. 5. As was done in Ref. 5, it was assumed 
that deflection would not affect the interaction equation provided that the ef- 
fect of deflection on the critical moment for the beam loaded by bending only 
is taken into account. 

The values of M', the critical moment without axial load, were first calcu- 
lated from the formula:(7 


M'= 


using L as the effective length in the weak direction, and where a2 = B yh2/4C 
and h = the distance between the centroids of the flanges. This gave values of 
M' = 17,840 in.-lb for the specimens tested with knife edges dry and M' = 
17,660 in.-lb for those tested with lubricated knife edges. These values do not 
include the effect of deflection. Since C is negligibly small compared to B, 
and Bo, the correction factor for deflection that appears in Eq. (1) reduces to: 


= = 1.160 
Bo lo 
By taking this correction into account, the critical moments become 20,680 
in.-lb and 20,480 in.-lb for the two cases. Fig. 5 shows the test results plot- 


ted in interaction form, both neglecting and including the effect of deflection 
on the value of M'. The curve marked “Timoshenko” represents the equation: 


(My*, 
t(yr) =! 


77\/ B,C a? (5) 


(6) 


which was derived in different form by Timoshenko for lateral buckling of 
deep rectangular beams, under simultaneous end load and uniform bending, 
and which was derived for I-section members by Johnston.(8) The curve 
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labeled “Goodier” represents the equation: 


p' ( ) 

where p is the polar radius of gyration of the section. This is the special 

case of a doubly symmetrical eccentrically loaded column of the more general 

solution by Goodier.(9,10) The remaining curve represents the equation: 


which is the same as Eq. (7), except for the factor (1 - a), which is included 
to account approximately for the increase in eccentricity of the load due to 
deflection in the plane of the web. 6 

It can be seen from an inspection of Fig. 5 that the test results closely fit 
Eq. (8) when the stiffening effect of deflection is included in the calculation of 
M'. Table I gives the results of the tests and also theoretical values of the 
ratio P/P"' calculated from Eq. (8) for each experimental eccentricity. Per- 
centage errors are given for the results both neglecting deflection and includ- 
ing its effects. The maximum difference between calculated and experimental 
values is 7.5 per cent when deflection is neglected, but only 2.3 per cent when 
the effect of deflection is included. 

Ref. 6 gave results of eccentrically loaded column tests on specimens simi- 
lar to those in the present investigation except that the column ends were fixed 
in the lateral direction. The differences between measured and computed 
critical loads in Ref. 6 are all within 3.5 per cent, and are all in the direction 
such that the correction factor for deflection which was developed for members 
with ends simply supported in the lateral direction would increase the dis- 
crepancy. Hence, it appears that this correction factor does not apply to 
members with laterally fixed ends. On the basis of these test results and the 
results of the theoretical study of rectangular beams, it appears that the ef- 
fect of deflection on the buckling strength of beams and beam-columns with 
ends laterally restrained is negligible. 

Since in practice the end conditions of beams and columns fall somewhere 
between the limits of simple support and full fixity, the effect of deflection on 
the critical moment is small, and its value in any particular case would be 
difficult to predict. It is, therefore, recommended that the effect of deflection 


on the critical moment not be included in the preparation of design rules for 
beams and beam-columns. 


CONC LUSIONS 


The following conclusions have been reached as a result of this study of the 
effect of deflection in the plane of the applied bending moment on the lateral 
buckling strength of beams and beam-columns: 


1. For beams of narrow, rectangular cross section subjected to uniform 
bending moments, the lateral buckling strength, including the effect of de- 
flection, is given by Eq. (1) when the ends are simply supported in the lateral 
direction and by Eq. (2) when the ends are fixed against lateral rotation. 

2. Beams of narrow, rectangular cross section may be appreciably 
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strengthened by the effect of deflection when the ends are simply supported in 
the lateral direction, but the deflection does not significantly affect the lateral 
buckling strength when the ends are fixed against rotation in the lateral di- 
rection. 

3. Experimental values of critical load for eccentrically loaded H-section 
columns with ends simply supported in the lateral direction show good 
agreement with calculated critical loads based on the interaction formula, Eq. 
(8), the value of critical bending moment in the interaction formula being de- 
termined from Eq. (23). The former equation takes account of the increase in 
bending moment resulting from the deflection of the member, and the latter 
equation includes the effect of the deflection on lateral buckling strength. 

4. Previous investigations have shown good agreement between measured 
and calculated critical loads for I-beams and for I- and H-section beam- 
columns(4,6) when 1), the ends were fixed in the lateral direction and 2) no ef- 
fect of deflection was considered except for the increase in bending moment 
due to deflection in the case of the beam columns. This result is in agreement 
with the conclusion, based on the analytical study of narrow, rectangular 
beams, that deflection has a negligible effect on the critical moment for later- 
al buckling when the ends are fixed in the lateral direction. 

5. Since in practice the end conditions of beams and columns fall some- 
where between the limits of simple support and full fixity, the effect of the 
deflection on the critical moment for lateral buckling is small, and its value 
in any particular case would be difficult to predict. It is, therefore, recom- 
mended that the effect of deflection not be included in the preparation of design 
rules for beams and that, in the case of beam-columns, only the increase in 


bending moment due to deflection be included, as previously recommended in 
Ref. 6. 


REFERENCES 


1. Michell, A. G. M., “Elastic Stability of Long Beams Under Transverse 
Forces,” The London, Edinburgh, and Dublin Philosophical Magazine and 
Journal of Science, Vol. XLVII, 5th Series, 1899, p. 298. 


2. Pettersson, Ove, “Combined Bending and Torsion of Beams of Monosym- 


metrical Cross Section,” Bulletin No. 10, Division of Technology, Stock- 
holm, Sweden, 1952, p. 55. 


3. Kerensky, O. A., Flint, A. R., and Brown, W. C., “The Basis for Design of 
Beams and Plate Girders in the Revised British Standard 153,” Proceed- 
ings, Institution of Civil Engineers, Part III, Vol. 5, 1956, p. 396. 


4. Hill, H. N. and Clark, J. W., “Lateral Buckling of Eccentrically Loaded I- 
Section Columns,” Transactions, ASCE, Vol. 116, 1951, p. 1179. 


5. Karol, Jacob, Discussion on Ref. 4, Transactions, ASCE, Vol. 116, 1951, 
p. 1192. 


6. Hill, H. N. and Clark, J. W., “Lateral Buckling of Eccentrically Loaded I- 


and H-Section Columns,” Proceedings of the First National Congress of 
Applied Mechanics, p. 407. 


7. Timoshenko, S., “Theory of Elastic Stability,” McGraw-Hill Book Company, 
Inc., New York, 1936, p. 261. 


1 
=. 
its, 
a 
% 
7 
5 
= 


ASCE EFFECT OF DEFLECTION 1596-9 


8. Johnston, Bruce, “Lateral Buckling of I-Section Column With Eccentric 


End Loads in Plane of the Web,” Journal of Applied Mechanics, December, 
1941, p. A-176. 


9. Goodier, J. N., “The Buckling of Compressed Bars by Torsion and 


Flexure,” Bulletin No. 27, Cornell University Experiment Station, Ithaca, 
New York, 1941. 


10. Goodier, J. N., “Flexural-Torsional Buckling of Bars of Open Section,” 


Bulletin No. 28, Cornell University Experiment Station, Ithaca, New York, 
1942. 


APPENDIX 


To determine the effect of deflection on the lateral buckling strength of a 
deep rectangular beam with ends fixed against lateral rotation, consider the 
equilibrium of any section m-m: 


The curvature 1/r is caused by deflection of the beam in the plane of the ap- 
plied moment M and is equal to M/Bz. The beam ends are fixed against ro- 

tation except in the plane of the moments. Considerations of equilibrium and 
symmetry of the beam require that the torques at the ends of the beam about 
the axis coinciding with the initial position of the longitudinal axis have zero 


magnitude. The equations of equilibrium of section m-m in the buckled con- 
figuration may be written: 
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AM-M,cos ($-#)=B, (F-4 (9) 


dy (10) 
sin (& - 

The right sides of Eqs. (9) and its represent the xn of St. Venant of 

circular beams bent out of their plane of curvature. 


Substituting M/Bz = 1/r, and eliminating v, one obtains the differential 
equation: 


Bt )(\- 
ds@ B,C (11) 
B M 
M(L- 
Letting 


By, C 2 
ac B, ge) =k 
the solution of Eq. (11) has the form: 


A=Asin ks +B cos ks Mi cos (5 - -s) (12) 
2 


When s = 0, Eq. (10) reduces to: 


M 
ds C Bo 2 (13) 


Evaluating A and B from the boundary conditions that # = 0 at the ends of the 
beam, substituting in Eq. (12), differentiating and combining with Eq. (13) 
gives: 


-k tan 


)tonM 


Be 2 


C Be (14) 
Since (M/Bz) (L/2) is small, replace its tangent with the angle. The angle 
kL/2 has a value very near 7 and, therefore tan kL/2 can be replaced by 
(kL/2 - 7) in Eq. (14). Solving the resulting equation for M, the expression 
for critical moment, M', is obtained: 


*The term (1/r dv/ds) was not included by Michell in his work.’ 
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M'= (15) 


The critical moment for an I-beam with laterally hinged ends subjected to 
uniform bending taking the effects of deflection into account can be similarly 
obtained. The equations of equilibrium are changed from Eqs. (9) and (10); 

M, is zero because of no end lateral restraint and the term B,h2 d3f is intro- 


4 
ie, duced to account for the resistance of the beam flanges to warping. 
i 
AM=B, (F - (16) 
af 
ds r ds 4 ds° .. 
Substituting M/B2 = 1/r and Pe aol V, one obtains the differential —_— 
equation: 
ds* ds Bo 
4c 
(18) 
and (1- 2) (1- 4, 
Bo’ d 
The solution of Eq. (18) is of the form: 
where m= 
and n= 
with the boundary conditions 
at s=O and s=L 
=O at s=O ands=L | 


ae 
q 
8 
q 
if 
4 


1596-12 EM 2 April, 1958 


The evaluation of the constants of integration from the boundary conditions 
shows the C2 = C3 = C4 = 0, so that 


B =C sin mz 


and sin mL=O 
The smallest root of this equation is 


(19) 


Substituting the values of m,—}-and i and solving for M, one obtains the ex- 
a 


pression for critical moment, 


The factor | , representing the effect of deflection 


B 
Be B 
for this case, is the same as the correction factor for the case of rectangular 
beam with the same end conditions. 
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CORRECTION FACTOR 


DEPTH/WIDTH RATIO 


Fig. | 


Effect of Deflection on Critical Moment of a 
Rectangular Beam Under Uniform Bending 
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AVERAGE DIMENSIONS AND CROSS-SECTION PROPERTIES OF SPECIMENS 


AVERAGE MECHANICAL PROPERTIES OF SPECIMENS 


Location Tensile Tensile Yield Cempressive Elongation 
Strength, Strength Yield Strength per cent 
psi 0.2 per cent 0.2 per cent in 2 in. 
offset, offset, 
psi psi 
me Flange 73 630 67 470 67 670 9.8 
' ig Web 74, 630 67 970 68 770 10.5 
Average 74, 130 67 720 68 220 10.1 


7 Tensile properties were determined from standard tensile specimens. 
(See "Standard Methods of Tension Testing of Metallic Materials," (E8-46), 
1949 Book of ASTM Standards). 


Compressive yield strengths were determined from 5/8-in. wide 
specimens of full thickness, (See "Tentative Method of Compression Testing 
of Metallic Materials in Sheet Form", (E78-4,9T), 1949 Book of ASTM Standards). 


Fig. 2 = AVERAGE DIMENSIONS AND PROPERTIES OF SPECIMENS 
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Determination of Eccentricity for 55-in. 22 
H-Section Columns with Knife Edge Ends 
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Fig. 5 


Interaction Relations for 55-in. 2-3/8" 
H-Section Columns with Knife-Edge Ends 
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INCREMENTAL COMPRESSION TEST FOR CEMENT RESEARCH#@ 


A. Hrennikoff,! M. ASCE 
(Proc. Paper 1604) 


This paper presents a description of a new compression test of cement and 
concrete which is so conducted that it separates the immediate deformation 
from creep and by that permits the determination of a number of elastic and 
inelastic characteristics describing in quantitative manner the stress-strain 


properties of material and throwing light on the internal mechanism of defor- 
mation and failure under load. 


In the course of experimental studies of the deformation and strength of 
cement and mortar the author has devised a special way of carrying out the 
compression test which results in a considerable amount of additional infor- 
mation not ordinarily obtainable. The compression test so conducted promis- 
es to be an important new research tool for investigation of mechanical 
properties of cement, mortar and concrete. Under the name of incremental 
compression test, it is described below together with the instrumentation and 
procedures used in connection with it. 

The stress-strain diagram of neat cement and concrete in compression is 
always presented as a curved line resembling a parabola with a steep initial 
part flattening out gradually toward the point of ultimate strength. A question 
naturally arises whether this non-linearity of the diagram is an indication of 
an early disruption of cement or a manifestation of some other cause. The 
new test partially answers this question. 

It is well-known that cement and concrete, as many other engineering ma- 
terials, are subject to creep, the word creep signifying a deformation requir- 
ing some time to develop. When material is loaded it becomes deformed 
simultaneously with the loading. This is the usual deformation which may be 


Note: Discussion open until September 1, 1958. To extend the closing date one month, 
a written request must be filed with the Executive Secretary, ASCE. Paper 1604 is 
part of the copyrighted Journal of the Engineering Mechanics Division, Proceedings 
of the American Society of Civil Engineers, Vol. 84, No. EM 2, April, 1958. 

Paper will be presented at a meeting of the ASCE, Portland, Ore., June, 
1958. 


Prof. of Civ. Eng., Univ. of British Columbia, Vancouver, Canada. 
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either elastic or inelastic. If however the load is kept constant the defor- 
mation normally continues to increase with time and this additional defor- 
mation is the creep. The new test allows the two kinds of deformation to be 
separated. 

The incremental compression test consists of loading the specimen very 
quickly in brief increments over parts of the compression range and in keep- 
ing the load constant after each rise for a long period during which the speci- 
men creeps. When there is no intention to carry the test to failure, the speci- 
men is unloaded in similar steps, each step followed by a period of back creep. 
The stress-strain curve of the test, consisting of the inclined loading and un- 
loading parts and horizontal ledges of creep, is plotted by the autographic 
stress-strain recorder. A special apparatus is also used for measuring later- 
al deformation. The procedure was first tried in the summer of 1954 and 
since then more than a hundred specimens have been tested in this manner. 

So far the tests have been performed only on standard cement and mortar 
cylinders 2" in diameter and 4" long. The 2" cubes are unsuitable for these 
tests for two reasons; firstly, because they have too little room for attachment 
of the strain recording apparatus; and secondly, because the stress condition 
in the cube is highly complicated and irregular. 

Ignoring for the time being the nonhomogeneity of the mortar or the cement 
material it is of course desirable for the stress condition in the specimen 
tested to be uniform, unidirectional compression. This ideal however is un- 
attainable, partly on account of imperfections of the ends and partly for the 
reason of friction on the end surfaces resulting from the difference in lateral 
deformations of the specimen and the contact surfaces of the testing machine. 
This effect is particularly pronounced in cubes where it creates a highly 
complicated three dimensional stress condition. 

Cylindrical specimens are also not perfect and in one respect are even in- 
ferior to the cubes; they must be capped on the upper end (with cement) for 
smoothness. This additional operation must be done very carefully. Further- 
more, the binding of this cap with the body of the specimen is never fully 
monolithic, Still, axial symmetry of cylinders and the more regular stress 
conditions resulting from it make them preferable to the cubes. For this 
reason they have been used in this investigation almost exclusively. 

Even well made cylinders are not free from the end effects in the course 
of testing and so a special manner of holding them is devised as shown in Fig. 
1. This arrangement involves steel caps 3" in diameter, provided with shallow 
recesses 2-1/16" in diameter, filled with hard rubber pads 1/8" thick, fitting 
tightly inside the recesses. The rubber smoothes the end irregularities and 
eliminates, at least partially, the end friction by deforming laterally in re- 
sponse to the lateral deformation of the cylinder, yet it is so hard (the index 
number 100 on the Shore rubber hardness meter) that its flow into the clear- 
ance recesses around the specimen is negligible even under compression 
stress as high as 15-20 kips/sq. in. Further improvement is effected by pour- 
ing liquid paraffin on the rubber disc prior to placement of the cylinder into 
the recess. The thin layer of paraffin formed provides lubrication between 
the rubber and the cylinder and fills minor irregularities. 

That the caps so designed are beneficial in reducing the undesirable end 
effects is manifested by the fact that in the majority of tests carried to failure 
the ends of the specimens remain intact until after the appearance in the 
stress-strain characteristics of signs of impending failure (to be mentioned 
later) thus indicating that the failure is not brought about by the end conditions. 
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In cases when the end defects are responsible for an early failure the latter 
is usually brought about by chipped edges (Fig. 2) causing the flow of rubber 
into the enlarged edge clearance space, thus dragging the adjacent cement 
portion outward by friction and cracking it off by tension on an inclined plane. 
The presence of paraffin lubrication is evidently useful in this situation. 

If a compression cylinder is supported by open rubber pads, unconfined in 
the cap recesses and devoid of paraffin lubrication, the friction forces de- 
veloped on the end of the cylinder by the outward flow of rubber are so great 
that the cylinder often splits in two along a diametral plane (Fig. 3) under a 
load less than a half of the usual. These manners of failure in the presence 
of rubber pads are different from the usual type of two-cone or one cone and 
split end failures observed in the cylinders with freely resting unprotected 
ends. 

The testing machine used in an incremental compression test must be of 
low inertia, and must be capable of quick loading and of maintaining the load 
at a constant level. A Baldwin, Tate, Emery universal testing machine of 
60,000 lbs. capacity has been found well suited for this purpose. The load 
increments of 5 and 10 kips may be attained in it in 5 to 10 seconds. A 
Baldwin-Emery load maintainer, although not indispensable and not fully ef- 
fective, helps to keep the load constant during the periods of creep. Skill is 
needed in starting and stopping the machine at the beginning and the end of 
each load increment but it is easily acquired by practice. 

The stress-strain graph, or rather the load-strain graph, is plotted by the 
high-magnification Baldwin stress-strain recorder, microformer type, with 
special compressometer attachment. The scale of longitudinal strain in the 
record may be as large as one inch to 500 micro-inches per inch. 

The diagram of the compressometer attachment is shown in Fig. 4 and the 
attachment in the photographs of Figs. 5, 6 and 7. The apparatus consists of 
a stand made of the side plates A and the base B (Fig. 4a). Two light hori- 
zontal frames C and Cy, provided with counterweights D, are hinged to the 
stand at the points E, The lower frames Cy carries the microformer assembly 
consisting of the coil F and the core G supported on the spring H, while the 
upper frame C, through the medium of the pin I, touches and depresses the 
core G. The right-hand sides of the frames C and C, make contact with the 
test specimen J by means of the pins K. For this purpose two horizontal rings 
L (Fig. 4b), slightly larger in diameter than the specimen, are hinged to it by 
the set screws M, with one pair of screws per ring. The rings carry vertical 
anvils N on which the pins K rest. The distance between the set screws M is 
the gauge length of the specimen. As the specimen shortens in the course of 
testing the rings L approach each other and so do the right-hand ends of the 
frames C and C, while their left-hand sides move apart. This induces the 
core G to emerge from the coil and thereby actuate the microformer and cause 
rotation of the recorder drum. No microformer action results when both 
rings L move down equal amounts. 

A special apparatus has been designed for measurement of lateral strains. 
As shown in Fig. 8, it consists of three solid segments A jointed by three bent 
iso-elastic strips B, each strip being fastened to the segments by small 
screws. The ring so formed is clamped on the cylinder by three large screws 
C passing through the mid-points of the segments. The iso-elastic springs B 
form the sensitive elements which bend as the specimen increases in diameter 
when it is subjected to axial compression. The bending strains in the iso- 
elastic springs are registered by short S.R. 4 strain gauges glued on the 
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inside and outside of all springs. All inside gauges and, separately, all out- 
side gauges are connected in series and are attached to the active and compen- 
sating terminals of the strain indicator whose readings are closely proportion- 
al to the lateral strains in the cylinder. The sensitivity of the apparatus is as 
high as one micro-inch per inch. The segments and screws are made of invar 
in order to reduce the error due to possible temperature variation. 

In the photograph of Fig. 7 the apparatus is shown mounted on the specimen 
together with the compressometer attachment. 

On the whole the apparatus has proved quite successful. Its weakness lies 
in the possibility of overstraining the iso-elastic strips by a sudden failure of 
the specimen. It is however expected that this defect will be eliminated by a 
modification in design. The apparatus can be easily adapted to any size of the 
cylindrical specimen. 

The significance of the incremental compression test becomes apparent 
from consideration of the stress-strain records and of the lateral deformation 
data of some of the tests. 

Figs. 9 and 10 represent two such records referring to the same cylinder 
of neat cement, 2" in diameter, which has been tested several times in suc- 
cession without bringing it to failure. The composition of the cylinder is 
0.003 : 0.24 : 1 by weight of pozzolith (a dispersing agent) water and cement, 
type 1, manufactured by the British Columbia Cement Company. The speci- 
men, 11 months’ old at the time of testing, had been kept continuously under 
spray in the moist closet and was tested moist. The specimen is subjected to 
3 incremental tests with a few days’ intervals apart. All these tests are 
carried to the load of 35 kips, or the stress of 11.15 kips/sq. in., and are con- 
ducted to the same time schedule. The times, in minutes, elapsed from the 
beginning of each test to different stages of it are stated along the graphs. 
Figs. 9 and 10 represent tests #1 and #3 respectively. 

As stated earlier, the stress-strain record consists of the inclined parts 
over which the load is raised or lowered very quickly and of the horizontal 
parts where the load is kept constant and the specimen is undergoing creep 
sometimes for a period of more than an hour. The most striking feature of 
the graphs is the comparative straightness and parallelism of the inclined 
parts indicating a nearly perfect proportionality between stress and strain and 
excellent elasticity of material, both in the course of loading and unloading. 

The time rates of creep, and the total amounts of creep on different hori- 
zontal steps, follow an easily recognizable regularity. During the loading part 
of the cycle they increase with the rise of the load level of creep, and during 
the unloading stage they increase as the load level returns to zero which 
means that the line depicting unloading appears similar to the line depicting 
loading rotated through 180°, 

The creep curves corresponding to several horizontal steps, when plotted 
on the basis of time, all present an appearance similar to Figs. 11 and 12. 
The time rates of creep in these curves are represented by the slopes. The 
difference between different curves is mostly quantitative. Fig. 11 refers to 
creep on the upper levels during the loading stage. The curve is steep at the 
origin, gets flatter with time, but is still fairly steep just before the load is 
raised further to a new level. A high rate of creep at the end of the curve may 
signify an impending failure. The lower levels in the course of loading and 
all levels during unloading have their creep curves similar to Fig. 12 charac- 
terized by a faster flattening and an asymptotic approach to a horizontal 
tangent. The creep rate at the origin of the curve on the zero load step, 
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following complete unloading, is very high, even higher than at the origin on 
the uppermost loading step. Some residual deformation remains after com- 
pletion of the test. 

Quantitatively the creep curves on all steps depend also on the time allowed 
for creep on the immediately preceding steps. For this reason a definite time 
schedule has been devised and adhered to in testing similar specimens for 
fair comparison of their creep characteristics. The schedule is somewhat 
arbitrary except for the requirement that the creep on each step continues un- 
til its rate is reduced to at least one tenth of the initial rate. 

It is interesting to note that in the records of Fig. 9 and 10, the immediate 
and creep deformations are separated almost completely. This may be demon- 
strated by examination of the creep rate just before and just after the last 
rise of the load from 30 to 35 kips in Fig. 9. Before this rise the specimen 
creeps about 40 micro-inches per inch in 10 minutes, or the creep rate is ap- 
proximately 4 micro-inches per inch per minute. At the top end of the rise 
the specimen creeps about 50 micro-inches per inch in the first minute and 
about 40 in the second minute. As the load is being raised the creep rate 
should be close to 4 - 5 micro-inches per inch per minute at the bottom and 
about 50 near the top with an average 1/2 (4 + 50) = 27 micro-inches per inch 
per minute. Since the rise is effected in 5 - 10 seconds, the amount of creep 
which becomes combined with the elastic deformation should be of the order 
(27)(10)/60 = 4.5m-in./in., which is almost unnoticeable on the scale of the 
record, 1 inch to 500™~-iN. /jn, This explains the straightness of the inclined 
parts of the graph. Some deviation from linearity may however be observed 
at the approach to zero load in the course of unloading, where the creep con- 
ditions are less favorable. Selection of greater load intervals between the 
periods of rest may result in a substantial intermixture of the elastic and 
creep deformations. 

When a specimen is tested again the second time following an identical 
time schedule, the new stress-strain record, being on the whole of the same 
type, displays some significant differences as well as similarities when com- 
pared to the first record. In further successive tests these features become 
more accentuated. This is demonstrated by comparison of test No. 3 (Fig. 10) 
with test No. 1 in Fig. 9. 

On repeated loading the new modulus of elasticity, both during loading and 
unloading, becomes somewhat smaller the average values in Fig. 10 being 
3960 kips/sq. in. on the way up and 3890 kips/sq. in. on the way down, com- 
pared to the values 4170 kips/sq. in. and 4150 kips/sq. in. in Fig. 9. 

The amounts and rates of creep during loading on all levels except the 
upper two are respectively mQch greater in Fig. 10 than in Fig. 9. The oppo- 
site is true with regard to the upper level of creep which is smaller in Fig. 

10 than in Fig. 9. The second from the top creep levels are approximately 
equal in both these graphs. The amounts and rates of creep on the way down 
are nearly identical in both Figs. 10 and 9 on all levels except the two bottom 
ones where they are longer in Fig. 10 than in Fig. 9. 

Even a casual inspection of the graphs reveals numerous minor irregu- 
larities in the shape of curves. Many of these are caused by faults in oper- 
ation and the imperfections of the apparatus such as: 


(1) The oscillations of the load maintainer. 
(2) Free motion of the recorder pen, corresponding to some play between 


the gear teeth forming part of the mechanical transmission of load 
indication. 
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(3) Rocking of the iron core inside the coil of the microformer resulting in 
error in strains. 

(4) Inertia of the moving parts associated with the strain indication. 

(5) The two loops observed in Fig. 10 in the vicinity of loads 30 and 35 kips 
are produced by the misuse of the controls. 


However all imperfections of the graphs cannot be explained in this manner. 
Some of these are undoubtedly caused by breakage and damage inside the 
specimen too small and too highly localized to lead at that stage to an overall 
failure. 

Some significant results are also obtained with reference to lateral defor- 
mations. These may be most conveniently described in terms of the variation 
of the coefficient of lateral deformation u defined as follows: 


4 f (lateral ) 


where the denominator and the numerator represent the corresponding in- 
crements of the longitudinal and lateral unit strains. 

The coefficient u is thus similar to Poisson’s ratio, except that it refers 
to particular increments of unit strains, corresponding to some definite inter- 
vals of time or increments of load, rather than to the full values of unit 
strains measured from the zero load condition. 

Fig. 13 represents the values of the coefficient “ plotted against the stress- 
strain record of Fig. 9. The graph shows that on the elastic parts uw is nearly 
constant and close to 0.2 both in the course of loading and unloading. Only a 
minor difference from this value is observed in the several consecutive tests. 

An entirely different situation exists with regard to the value of uw under the 
conditions of creep. At the lower levels of creep on the way up the value of u 
is close to zero and as the load level rises it increases to the value approach- 
ing and even surpassing 0.2. At each particular level of creep the quantity yu 
tends to remain constant. Some scattering of this quantity is however ob- 
served from one time interval to another on the same horizontal steps in view 
of the difficulty of measuring accurately the small strains involved. 

In the course of unloading the values of uw associated with creep are more 
constant than in the course of loading and are intermediate between zero and 
0.2, as is evidenced by Fig. 13. 

Failure of a moist cement specimen, occurring usually during creep, is 
apparently related to strain rather than to stress, and its approach is indicat- 
ed by several signs: 


(1) The longitudinal strain at failure is usually greater than 5000 micro- 
inches per inch and sometimes it even exceeds 7000. 

(2) The time rate of longitudinal creep just before failure is high and re- 
mains constant or even tends to increase on the same creep level. This 
is in contrast with the behaviour at the lower levels where the creep 
rate gradually decreases. 

(3) The value of during creep at the approach of failure is more than 
0.25 and it tends to increase at the same creep level. 


The behaviour of cement under compression described so far is exhibited 
by the specimens which are moist at the time of testing. A very different 
stress-strain record is produced by the dry specimens. This is illustrated 
in the graph of Fig. 14 referring to a specimen of the same composition as the 
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moist cylinder considered above. This cylinder was dried before testing to a 
constant weight in an oven at 200° F. In the course of drying the specimen lost 
slightly more than 50% of the weight of water used in its making. 

The incremental stress-strain graph of a dry specimen consists mostly of 
straight inclined lines corresponding to quick loading and unloading. The hori- 
zontal steps of creep between these rises and falls are very short and nearly 
non-existent at the lower levels of loading and the upper levels of unloading. 
In view of the virtual absence of creep the time schedule of the test has been 
much shortened compared to that of the moist specimens. Unlike the moist 
specimen, it is sometimes difficult to see any regularity in the relative magni- 
tudes of the creep steps at different levels, and on the whole the record is 
less regular than that of a comparable moist specimen. The modulus of 
elasticity, also nearly constant over the whole range of loading, is nearly the 
same as in a similar moist specimen. This is also true for the value of 
over the elastic part of the graph, which is again close to 0.2. Over the creep 
parts of the graph the values of uw are difficult to determine accurately in view 
of the smallness of deformations but they appear to be close to the values ob- 
served in the moist specimens. On repeated loading of a dry specimen, as is 
illustrated in Fig. 15, creep becomes even smaller but the modulus of elastici- 
ty remains unchanged. 

Both in Fig. 14 and 15 the stress-strain graph has a peculiar characteristic 
feature near the zero load. The graph starts off at the beginning of the test 
and ends at its conclusion with some comparatively flat portions, and further- 
more, when the graph becomes steep and straight in the course of the first 
rise and fall the corresponding moduli of elasticity are appreciably smaller 
than in the subsequent parts. The reason for this anamaly is the presence in 
the cylinder of fine transverse cracks which form in the course of drying and 
are virtually unavoidable. Closure of these cracks at the commencement of 
loading and their opening and unloading is responsible for the additional strain 
deformation present. 

The mechanical behaviour of moist and dry cement described here is quite 
general, and it has been observed on many dozens of cement specimens of 
different kinds. Age, water-cement ratio and coarseness of grains produce 
some significant modifications in the stress-strain records, but the general 
type of action remains unchanged. 

The preceding discussion shows how the incremental compression test 
separates creep from the immediate deformation and allows the numerical 
values of various parameters describing different aspects of stress-strain 
behaviour of cement to be determined. The several parameters and charac- 
teristics so found may be enumerated as follows: 


(1) The modulus of elasticity and its constancy (or lack of it) throughout the 
test. 


(2) The law of variation of longitudinal deformation with time during creep 
at each load level and from level to level. 
(3) The coefficient of lateral deformation u 


(a) during elastic deformation 
(b) during creep at different levels of load. 
(4) The above quantities may be found both in the course of loading and in 
the course of unloading. 
(5) The same quantities may be found during several consecutive tests on 
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the same specimen. The differences reflect the effects of previous 
loadings. 

(6) Significant similarities and differences in the stress-strain characteris- 
tics are observed on comparison of the moist and dry specimens of the 
same composition and of compositions different with regard to the 
water-cement ratio and the fineness or coarseness of grains of cement. 

(7) The signs of approaching failure in terms of stress-strain characteris- 
tics may be recognized and studied. 


The applicability of the incremental test is not limited to neat cement but may 
be extended to mortar and concrete. 

Several aspects of the mechanical behaviour of cement presented here are 
the manifestations of the physical and chemical factors at work which pre- 
sumably may, with the assistance of imagination and good luck, be brought to 
light by deductive reasoning. Thus, the incremental test provides a new 
avenue of approach in the study of cement and concrete. What actions and what 
constituents in cement paste are responsible for its stress-strain behaviour 
is a fascinating problem. Some measure of progress towards its solution has 
been made in the course of these studies but the subject is too wide for an 
adequate discussion at this stage and is reserved for the future. 


The experimental work described here has been conducted in the Materials 
Testing Laboratory, Department of Civil Engineering, The University of 
British Columbia. The work has been given generous financial assistance by 


the Department of Building Research, The National Research Council of 
Canada, The author wishes to express his deep appreciation for this as- 
sistance to the President of the National Research Council, Dr. E. W. R. 
Steacie, to the Director of Building Research, Dr. R. F. Legget, and to the As- 
sistant Director of Building Research, Dr. N. B. Hutcheon. 

The author also wishes to thank Mr. J. F. Muir, Head, Department of Civil 
Engineering, The University of British Columbia for his cooperation in mak- 
ing available to the author the laboratory and shop facilities of the Department. 
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WIND INDUCED VIBRATION OF CYLINDRICAL STRUCTURES# 


Closure by Joseph Penzien 


JOSEPH PENZIEN,* J.M. ASCE.—The writer wishes to thank both Mr. Au 
and Mr. Roshko for their informative discussions on the vibration of cylin- 
drical structures subjected to wind loading. The increased interest in this 
field of research recently is indeed gratifying as the problem of stack vibra- 
tion is of growing concern to industry and it is perhaps safe to say the pro- 
fession does not know with certainty how to design a tall steel stack with 
proper consideration of the dynamic behavior under wind conditions. When 
one considers the high cost of such a structure, it is obvious additional re- 
search is desirable in this field. 

The test results presented by the writer show the critical range of Rey- 
nold’s numbers somewhat below that found by Relf and Summons.(1) As 
pointed out by Mr. Au, the results of Relf and Summons are very similar to 
those results found by Delany and Sorensen.(2) This difference would seem 
to indicate a higher level of turbulence in the wind tunnel used by the writer 
than in those used by the above mentioned investigators. 

The writer is especially grateful to Mr. Roshko for his discussion on the 
intermittent changes in flow which may produce “bursts” of vibration at 
natural frequency in the critical range. While these “bursts” of vibration are 
Similar to those observed on stacks well above the critical range, it is not 
known if the basic phenomenon causing such response in both ranges is the 
same. Further research should be concentrated on this basic aerodynamic 
behavior. 

Mr. Roshko’s criticism regarding lack of information on amplitude of re- 
sponse and magnitude of aerodynamic lift is favorably received. To correct 
this oversight, the spring constant “k” for the system used in each test is 
given in the subsequent table. However, the writer did not attempt to discuss 
the actual magnitude of aerodynamic lift and amplitude of vibration in his pa- 
per aS maximum response depends not only on the aerodynamic lift but also 
on the time duration of the resonance condition. The real question is under 
what conditions does the aerodynamic forcing function stay in resonance with 
the system sufficiently long to produce excessive vibrations. The results 
presented show the resonance condition to exist for all wind velocities in the 
critical range which is similar to that behavior observed on stacks above the 
critical range. This condition is unstable which results in the “bursts” of 
vibration. However, under certain conditions of structural stiffness, struc- 
tural damping, wind velocity, etc., the resonance condition becomes relative - 
ly stable; thus producing excessive vibrations which may result in structural 


a. Proc. Paper 1141, January, 1957, by Joseph Penzien. 
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failure. It is therefore the basic aerodynamic phenomenon associated with 
this stable condition which is of practical importance and needs further in- 
vestigation. 


Spring Constants k - lbs/in. 


Test 
No. k 


890 
396 
172 
113 
43.0 
955 
405 
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No. No. 
No. k No. k 

1 2830 8 15 149 ra 

2 1120 9 16 116 a 

3 238 10 17 41.7 a 

4 153 11 18 314 : 

5 86 12 19 126 : 

6 33.8 13 20 37.8 : 

7 1475 14 
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BENDING OF A RECTANGULAR PLATE WITH ONE FREE EDGE@ 


Closure by W. Nachbar 


W. NACHBAR.!—One purpose of the paper was to present a considerable 
quantity of accurate numerical data against which the predictions of approxi- 
mate, analytical or finite-difference methods for plates may be checked. 

Mr. Silverman has made excellent use of this in demonstrating that his 
method is a useful one for obtaining approximate, closed-form solutions to 
the problem under discussion. Although the original paper dealt only with 
the exact solution, it may be noted further that useful, approximate, closed- 
form expressions for the various quantities listed in Eqs. (26) through (36) 
can also be found by substituting into these equations the series given in Eqs. 
(24) and (25) and then neglecting the terms involving high powers of exp (-u). 


a. Proc. Paper 1196, April, 1957, by W. Nachbar. 
1. Research Branch, Missile Systems Div., Lockheed Aircraft Corp., Palo 
Alto, Calif. 
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DEFLECTIONS OF STRUCTURES IN THE INELASTIC RANGE@4 
Closure by Kurt H. Gerstle 


KURT H. GERSTLE, | J.M. ASCE.—Thanks are due to Mr. Percy for his 
valuable discussion, in which the method for obtaining deformation coeffi- 
cients is generalized to apply to any type of moment-curvature relationship. 
With the expressions presented by Mr. Percy, it is possible to state informa- 
tion regarding the behavior of I-section and other beams, to take into con- 
sideration the effects of strain hardening of mild steel, and to predict de- 
formations of members of aluminum and other materials without a definite 
yield range. 

The technique of interpreting the coefficients for beam segments under 
linear moment variation in terms of the properties of the moment curvature 
diagram is useful, it is however believed that the definitions of the quantities 
Pg and Qg got interchanged. 

The reference to available literature on methods of obtaining moment 


curvature relationships is appreciated. Opportunity is taken here to point 
out errata of the original paper: 


Page 2, Line 28: (T= Sypt [c* $) 


Page 17, Line 1: “The frame shown in Fig. 12 and 13 is indeterminate to 


the second degree. Two moments must be estimated 
Simultaneously, such that...” 


a. Proc. Paper 1290, July, 1957, by Kurt H. Gerstle. 
1. Instr., Univ. of Colorado, Boulder, Colo. 


2 
ASCE 1619-7 ic, 
¥ 
2 
“has 
ofa 
4 
= 
a 
q 
. 
Me 
4 a 


‘ 
ag? 
al 


ASCE 1619-9 


DEMONSTRATIONS OF PLASTIC BEHAVIOUR OF STEEL FRAMES4 
Discussion by Zdenek Sobotka 


ZDENEK SOBOTKA.!—The authors have accomplished an extremely 
valuable job of demonstrating the main characteristics of the plastic behavior 
of steel frames and of confirming the principles of limit design. Particular- 
ly, the presented results of full scale tests are of great importance for limit 
analysis of steel structures. The authors’ tests prove in accord with former 
experimental and theoretical investigations the independence of plastic load- 
carrying capacity on the distribution of internal stresses in the critical sec- 
tions where the plastic hinges develop. This may be explained by the mutual 
equilibrium of these stresses. 

This discussion is intended to examine more in detail the effect of shear 
on reducing the plastic bending capacity as mentioned on page 1390-21. 

The derivation of the interaction curve and of interaction equation repre- 
senting all possible limit combinations of bending moment and shear force in 
I-beams may be based on the assumption that the flanges transfer only bend- 
ing and none shear, while in the web act uniformely distributed bending 
stress o and shear stress T, as shown on Fig. 1. This assumption concides 
with that expressed by the formula (1) on page 1390-21. 

The flow condition 


1 
, 


which, for the web, may be written in the form 


w 


1 e 1 2 t (2) 
leads after multiplying by the plastic section modulus of the web, 4 A. , to 
the following equation of the cricle 
2 
r)2 va 


a. Proc. Paper 1390, October, 1957, by H. M. Nelson, D. T. Wright and J. W. 
Dolphin. 

1. Docent of Strength of Materials and Theory of Elasticity, Technical Univ., 
Prague, Czechoslovakia. 
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in the coordinates of and V 


q 
In the preceding equations 
M is the plastic moment of web, 
~ 4 y 


M,* the plastic moment of web reduced due to shear, 


V the shear force. 
The notations of cross-section dimensions may be seen from Fig. 1. 
The interaction curve of the I-beam for the limiting state of plastic bend- 


ing combined with the shear is plotted in the coordinates of M and Sel? 
on Fig. 2. The straight line BC represents the load-carrying capacity of 
flanges and the arch CD that of the web. 


The point A represents the limit combination of the bending moment M, al 
and of the shear force V, which may be analytically expressed by 


Mp+M, sing , (4) 


te) Oy (5) 


there is full plastic moment of flanges. 


The value of sin @ may be obtained from geometrical relations of the in- 
teraction curve, as follows 


f 
= 6) 
( 


4 


= te Xe (7) 
After some transformations, we get from Eq. (6) the relation 


sing = 


1 
tga 1+ tg% of 


and, after substituting it in Eq. (4), the general formula for the plastic bend- 
ing moment of the I-beam reduced by the effect of shear 


Me 
= My | tea 1+ ten - (9) 


where 
sin — 
Ww Vv 3 
where 

— 

M, 

f 
Ww 
r 
M 
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One may examine, for example the reduced load-carrying capacity of the 
propped cantilever with single vertical load (Fig. 3). The full unreduced load 
at collapse is on page 1390-9 given by 


L+b 

W = M,, 

The shear forces in the two parts AC and CB are Wy and W9 respectively. 

Using the method of forces, we have from the limit plastic equilibrium of the 
two parts of the propped cantilever 


r 
(10) 


from which 


ri 2 2) (12) 
q = 
| 


* 


For bee 5 , it follows from Eqs. (10) and (11) that Ww, > Wo: Then, the 
reducing effect of shear in the part AC is greater than that in the part CB, 


while the reduced bending moment M., is in the point C the same for both 


parts of the beam as follows from statics. 
Introducing the shear force V = W, from Ey. (10) in Eq. (7) leads to 


(13) 


After substituting the value of tga@ from Eq. (13) in Eq. (12), we get the 
reduced ultimate load 


a 2 
3M 
{p+ 4a” + La 


From (9) it follows, that the reducing effect of shear increases with the 
decreasing value of tga , or with increasing value of clear web depth dy and 
with the decreasing length a between two plastic hinges, as may be seen from 
Eq. (13). 

Thus, in the case of short beams with great depth of the web, the effect of 
shear on plastic load-carrying capacity is comparatively large and must be 
taken into account. This was confirmed by some results of tests of beams 
with high webs. 
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PLASTIC DESIGN OF COVER PLATED CONTINUOUS BEAMS# 


DOUGLAS T. WRIGHT.!—Messrs. Popov and Willis, in the synopsis of 
their paper, say “(the) situation (in plastic design) is entirely different from 
that found in elastic design where, theoretically, once the type of cross- 
section is fixed there is only one solution of a given problem.” And, again, 
on the third page of the paper, the authors say that for a given loading con- 
figuration and a given variation in section, there would only be a single 
unique elastic solution, whereas, ‘for the same conditions an infinite number 
of designs are possible according to plastic theory.” Although by reading 
between the lines, it is possible to understand what the authors were trying 
to say, the statements quoted above, taken at face value, are simply not true: 

For a given fixed loading configuration, (all loads proportional to P) anda 
given configuration of strength (all moments proportional to Mp) there is only 
a Single unique load factor (A) which will produce collapse according to the 
plastic theory.(1) By elastic theory there is also only a single unique solu- 
tion once loadings and sections are defined. By varying sections in the de- 
sign operation, an infinite number of solutions are possible both with elastic 
and plastic theories. 

Thus in these respects, at least, elastic and plastic methods are in fact 
identical in nature contrary to the authors’ statements. The significant dif- 
ference between elastic and plastic theories is that by elastic theory the 
limiting capacity is reached when a single section reaches yield whereas by 
plastic theory the limiting (fully plastic) capacity is not reached until hinges 
form at (n + 1) sections. (n being the number of indeterminacies). 

A very important practical difference between elastic and plastic designs 
of sections with varying capacities (i.e. cover-plated not noted by the 
authors) is that trial and error methods must be used in elastic design whilst 
direct design is possible by the plastic method. 

Later in their paper Messrs. Popov and Willis suggest that simple plastic 
theory does not satisfactorily estimate the full strength of the beams they 
tested. The test loadings used provided discontinuities in the bending moment 
diagrams at all points of maximum moment, which in turn lead to short hinge 
lengths, high strains, and early strain hardening. With distributed loadings, 
(which are of course more common in practice) strain hardening is not so 
important, and the simple plastic theory provides better estimates of 
strength.(2) Methods are in fact available by which the effects of strain har- 
dening can be readily considered. (3) 


a. Proc. Paper 1495, January, 1958, by E. P. Popov and J. A. Willis. 
1. Asst. Prof. of Civ. Eng., Queen’s Univ., Kingston, Ontario, Canada. 
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